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ABSTRACT 


This thesis explores the physieal and mathematical limitations of two common attitude 
control systems: one based on reaction wheels and another based on control moment 
gyroscopes (CMGs). The dynamics are derived from first principles and control 
algorithms for achieving maximum reaction wheel potential are discussed. The shaped 
eigenaxis input is utilized to establish baseline maneuver performance. A time-optimal 
shaped input is introduced and implemented in a feedback setting, subject to the 
limitations of the Moore-Penrose pseudo-inverse control allocation. Finally, a feed¬ 
forward plus feedback controller is introduced to implement the time-optimal torque 
inputs directly to the reaction wheels. This obviates the need for the pseudo-inverse 
control allocation, and therefore exploits the total capacities of both the reaction wheel 
momentum envelope and torque envelope. These reaction wheel control approaches are 
compared with CMG performance to establish spacecraft size and slew parameters that 
make the use of reaction wheels a reasonable choice. 
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I. INTRODUCTION 


A. MOTIVATION 

Spacecraft in low earth orbit (LEO) traverse around the earth at very high rates. 
For instance, a spacecraft in a circular orbit at 200 nautical miles altitude travels at speeds 
in excess of 17,000 mph. This equates to approximately 4 degrees of the orbit traveled 
every minute, completing a full orbit in 90 minutes. From the point of view of an imaging 
satellite, the spacecraft may only be able to see a ground target for 3-4 minutes, but this 
amount of time requires the orbit to pass directly over the target. Because much 
resolution is lost when taking low angle photographs, only the 20 degrees off-nadir is 
conventionally used [1, Table I-l], [2]. The off-nadir constraint significantly reduces the 
amount of time the target is within the line of sight. Depending on how far left or right of 
the orbit track the target is, the satellite may only have seconds to obtain an image. The 
problem is further complicated when multiple targets he along that small target zone. 

This spaee-based problem is even passed forward to the military ground 
commander. The battlefield owner requires timely battlespace situational awareness for 
his mission but may not have the state-sponsored weight in order to gain access to 
national assets within a timely manner. The Army’s new imaging satellite concept. 
Kestrel Eye, is a new approach to fill this void. The intent is to put many small, 
inexpensive satellites into orbit to obtain images and transmit them to the battlespace 
owner more quickly. These satellites must be agile to get multiple images per target zone 
and responsive to last-minute request ehanges. Therefore, conventional slew rate 
limitation of about 2°/sec must be overcome [2]. The meehanisms that induee rotation 
(momentum exchange devices) must therefore be used to their highest potential. In this 
way, the spacecraft can move from target to target as quickly as possible to minimize 
otherwise unproduetive slew time. This thesis explores teehniques to maximize actuator 
performance and minimize slew time without any hardware changes. 
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B, BACKGROUND 

Control engineers typieally employ two different types of eleotro-meohanieal 
momentum exehange deviees (MEDs) to perform spaeeeraft attitude eontrol. The first is 
reaetion wheels, whieh eause a reaetionary body torque and momentum by eounter- 
rotating a small rotor. Reaetion wheels are simple to eontrol and have a high momentum 
eapaeity, whieh translates to high angular veloeity rates for the spaeeeraft. The seeond 
MED is a control moment gyroscope (CMG), which houses a spinning wheel within an 
intermediate gimbal frame. The gimbal frame is allowed to rotate, changing the 
orientation of the wheel’s momentum vector. The gimbal torque input is in turn amplified 
(as the product of the gimbal rate and the magnitude of the momentum vector). The high 
torque output of CMGs is complicated by mechanical complexity of multiple rotating 
frames and the presence of zero-output singularities that must be avoided [3]. The 
elevated torque, and therefore higher acceleration of the spacecraft body, has led CMGs 
to become the de facto solution to generate agile spacecraft. Although reaction wheels are 
more prevalent [4] (Appendix B), cost less, are mechanically simpler, weigh less 
(Appendix B), and are easier to control, the torque output of CMGs simply cannot be 
matched by reaction wheels. 

A spacecraft slew is the reorientation from one orientation in three dimensional 
space to another. Eor a basic slew, the initial and final spacecraft orientations are given in 
terms of their quaternion reference points. To reconcile the difference, the most basic 
approach is to treat the difference in quaternion vectors as a single arc. The arc of 
minimum distance rotates about a single axis, known as the Euler axis [5, p. 155] or 
eigenaxis [6, pp. 329-338]. A basic approach of control therefore is to treat the three- 
dimensional reorientation problem as a one-dimensional single angle slew along this arc, 
about the eigenaxis. The eigenaxis attitude control maneuver requires a single quaternion 
to be called for and error is regulated to zero as the spacecraft gets closer to the desired 
quaternion. This maneuver can be done as a step input, where the error is constantly 
corrected for by the controller [7]-[9]. The step input approach is simple to implement 
and stability is dependent upon only a few gains. Thus, the approach is popular in 
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spacecraft attitude control [10, pp. 212-267] [5, pp. 152-194] [6, pp. 403-452] [11, 
pp. 351-430]. 

One of the challenges of the step input approach is gain selection. The gains 
designed for large maneuvers may act slowly for small maneuvers, causing poor 
performance. Another approach is to use a shaped trajectory. Instead of passing a single 
final quaternion to the controller, a quaternion path for the maneuver is calculated and a 
set of discrete intermediate quaternions is passed as a target trajectory. Here, the gains are 
tuned to give acceptable performance and their values are agnostic to the slew size. 
Creamer et al. [12] developed a shaped bang-coast-bang maneuver about the eigenaxis, 
used on board the Clementine spacecraft. Time-optimal quaternion shaping was used on 
the TRACE spacecraft [13]. Additionally, [14] suggests shaping both trajectory and body 
rates and using these as inputs to the spacecraft attitude control (ACS) system. 

A typical ACS uses four or more wheels for redundancy. However, the use of four 
or more wheels presents a control allocation problem. This is caused when torque or 
momentum required in the body frame must be produced by the redundant set of MEDs. 
A common method of torque and momentum allocation is to use the Moore-Penrose 
pseudo-inverse, which provides a least squares (E2) solution [5, p. 169] [6, p. 440] [11, 
p. 416] [15]. According to this approach, the full magnitude of the requested body torque 
may not be available in every direction. The minimum norm solution is devised, 
providing a unique allocation of wheel torques [16, pp. 49-57]. The pseudo-inverse does 
not take physical wheel limitations into account and may request more torque than is 
available. Once any wheel is saturated, then the remaining allocated torque is simply lost. 
This is because to remain on eigenaxis, the integrity of the least squares solution must be 
maintained. Thereby, proportional saturation across all wheels is necessary. Because of 
the nature of the Moore Penrose pseudo-inverse, maximum momentum and torque will 
be unavailable about certain axes. The pseudo-inverse is described within [15] in great 
detail, providing detailed illustrations of the E2 solution in three-dimensional space given 
various reaction wheel array configurations. 

A way to bypass the E2 limitation is to operate the ACS in the open-loop 

(feedforward). One approach is to calculate the optimal torque solution and command the 
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wheels direetly. This is far from eonventional thinking, beeause the individual wheel 
torque solution is not generally available. The shape of the torque signal ean be derived 
over the eourse of the maneuver. This feedforward eoneept has been around for deeades 
[17]-[19], but the approaeh is susceptible to instability due to external disturbances and 
uncertainties of the spacecraft inertia or initial states [9]. For the solution to work 
properly, all states and parameters must be established exactly prior to reorientation, with 
no uncertainties throughout the maneuver. This ideal scenario is unlikely, so a merger of 
open-loop and closed-loop (feedback) must be used to improve stability and performance. 
Closed-loop implementation of open-loop maneuvers was flown on TRACE [20] with 
great success, but the performance gains were somewhat limited by the spacecraft 
preexisting software. Based on the literature available, it appears that industry may not be 
utilizing reaction wheels to their full potential, fostering a greater disparity between the 
capabilities of reaction wheels and CMGs. 

C. LIMITATIONS ON AGILITY 

In practice, a spacecraft slew rate is restricted for a variety of reasons. The most 
obvious limitations are due to the physical limitations of the MED array. Eor a reaction 
wheel system, an individual wheel is limited on how fast it can spin, thus the momentum 
it can produce. Positive rotation at maximum angular velocity results in a maximum 
momentum vector, while counter-rotating results in a momentum vector of equal 
magnitude and opposite direction. Given a skewed array of several reaction wheels, the 
vector sum of the maximum momentum vectors in any arbitrary direction creates the 
momentum envelope (Chapter IV). The shape of the reaction wheel momentum envelope 
is a polyhedron [15], involving sharp edges and vertices. This signifies that more 
momentum is available in certain directions than others. This momentum envelope 
produces the physical rate limit of the system about any given axis. 

Many attitude control systems involve restricting the axis of rotation to the 
eigenaxis, thereby minimizing the angle of rotation. Since the momentum space is not 
spherically uniform, some axes allow faster rotation rates. This causes additional 
problems with slew planning. Not only are two quaternions in space given, but there is 
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also an eigenaxis-specific slew rate. Not to mention the inerease in eigenaxis-speeifie 
slew acceleration (derived from the torque envelope, Chapter IV) to be concerned about. 

When operators want an imaging satellite to obtain an image, the rate and 
acceleration variations over 47i steradians require them to focus on the temporal aspect of 
the maneuver. Typical questions are; Given the slew angle and specific eigenaxis, will 
the spacecraft arrive at the desired quaternion at the correct time? Is there another axis 
the spacecraft can slew about more quickly to image something else? Even with 
automation, trajectory planning becomes a time-consuming task. It may be far easier 
from an operational point of view to operate within a single capability, even if it means a 
loss of performance. This procedural limitation further reduces the operable reaction 
wheel momentum space polyhedron to the largest inscribed sphere, which provides 
guaranteed capability achievable in every direction. A simple example is a sphere 
inscribed inside of a cube (Figure 1). The envelope (cube) restricts the magnitude of the 
vector to 1.0 in some directions, while the maximum magnitude is 1.73 at the vertices. 
This amounts to a non-uniform momentum space, offering values of momentum about 
some axes that are considerably less than others. If operating at the envelope, the motion 
is difficult to predict without correlating the specific eigenaxis with the available 
momentum. Therefore, the engineer might design to the lowest magnitude (sphere) and 
provide this momentum limit to the planners as their planning constraint for travel about 
any axis. The same logic can be applied for torque limitations, but the main point here is 
that some performance is lost due to the need to simplify other aspects of spacecraft 
operations. 
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Max Sphere = 1.0 
Max Envelope = 1.73 


Figure 1. Sphere Inseribed within a Cube 

D. OBJECTIVES AND SCOPE OF THIS THESIS 

This thesis is an attempt to elose some of the performance gaps between reaction 
wheels and CMGs, with the hope of identifying some key regions where reaction wheels 
may actually outperform CMGs. The true capabilities of reaction wheels are therefore 
explored to determine whether or not such an improvement is even feasible. The shaped 
trajectory used on the Clementine spacecraft [12] is developed (Chapter III) and analyzed 
(Chapter V). This maneuver becomes the basis for comparison of two different reaction 
wheel control algorithms that were developed as part of this thesis. 

The main development is the time-optimal maneuver, which merges shaped 
torque feedforward with shaped quaternion and spacecraft body rate feedback (Chapter 
V). This approach avoids the mathematical limitations of pseudo-inverse control 
allocation by applying torque signals directly to the reaction wheels. The caveat of using 
feedforward plus feedback is that a margin must be built in to account for uncertainties. 
This adds to the body of knowledge because the use of optimal quaternion and body rate 
trajectories has the ability to correct for spacecraft parameter uncertainties and 
disturbances. An approach for adding system margin based on a level of uncertainty is 
discussed. Optimal solutions are solved using DIDO [21] and then propagated through a 
computer simulation of the spacecraft dynamics to verily their feasibility. 
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The second reaction wheel enhancement is an extension of the time-optimal 
trajectory, but limited to the pseudo-inverse space. This method again breaks away from 
the eigenaxis, but operates entirety within the confines of the L2 volume. This allows the 
control engineer to upload the time-optimal quaternion trajectory instead of simply the 
eigenaxis quaternion trajectory [12]. Since the only change is the shape of the trajectory, 
it is not very intrusive on the current employment of a standard shaped quaternion 
trajectory implementation. This method will be introduced at the end of Chapter V. 

To aid in the setup of the problem, the spacecraft, reaction wheel, and CMG 
dynamics will be derived from first principles. This is done in Chapters II, IV, and VI, 
respectively. The equations of motion are then implemented within Simulink, which will 
act as a computer simulation of the actual nonlinear spacecraft dynamics. Various control 
algorithms will be applied to the same simulation to retain the structure of the model. 
Either the reaction wheel or CMG simulations can be selected as the current MED device 
within the model. Simulink also offers a graphical depiction of mathematical changes to 
the full dynamic setup, whereas changes can be hidden within loops if utilizing 
MATLAB alone. 

The sphere concept is applied to the time-optimal case. Eor a series of time- 
optimal maneuvers between any two quaternions, it is shown how the minimum 
guaranteed momentum and torque can be derived for any direction of travel about any 
effective eigenaxis. This adds the ability for the planner to devise a single slew rate over 
any rotation angle, because the sphere remains independent of rotation axis and angle 
size. The effective sphere increases in size twice from the smallest inscribed sphere, once 
for each of the two optimal control techniques mentioned above. Therefore, two levels of 
ACS improvement are given as part of this work. Chapter V presents the concept of 
effective eigenaxis for an optimal maneuver and the details of the derivation. 

A new technique for visualizing slews is introduced in Chapter VIE This method 
utilizes the normalized momentum and torque envelopes as the surface to display 
momentum and torque vectors throughout the maneuver. In the literature review, 
eigenaxis and quaternion traces were found, which only accounts for the path traveled. 

This new visualization tool allows the behavior of the rate and acceleration to be 
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observed with respect to the available momentum and torque capacity. The tool 
demonstrates fundamental features of the optimal control maneuvers. 

Finally, the performance of all three reaction wheel control approaches are 
compared against that of a like-sized CMG in Chapter VII. Data from the market research 
from over 75 reaction wheels and CMGs is used to generate trends and differences 
between these different systems. A simple method of comparing reaction wheels and 
CMGs is developed. Momentum and torque capacities are compared with slew angle 
sizes to determine regions where reaction wheels outperform CMGs under the various 
control algorithms tested. Concluding remarks and suggestions for future work are given 
in Chapter VIII. 
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II. SPACECRAFT KINEMATICS AND DYNAMICS 


The purpose of this ehapter is to develop the spaeeeraft kinematies and dynamies. 
Spaeeeraft attitude kinematies will be modeled using quaternions. Spaeeeraft dynamies 
will be derived from first principles, showing the dynamic relationship between the 
momentum exchange device (MED) action and the spacecraft body reaction. The MED 
assembly (subscript medA in this chapter) is any generic arrangement of momentum 
exchange devices. It can be 3, 4, or n reaction wheels, n CMGs, or any combination of 
thereof. 

A. MODELING A SPACECRAFT 

The spacecraft, MED assembly, and feedback can be modeled within Simulink, 
shown in Eigure 2. 


Orbital Motion 



Eigure 2. Spacecraft Simulation Model 


Given a desired (commanded) quaternion and a current quaternion state, a 
quaternion error vector is produced. The controller attempts to drive that error to zero by 
producing an appropriate torque on the spacecraft. The MED Assembly attempts to 
produce that torque. The Spacecraft Dynamics model the spacecraft response, with torque 
disturbances {Td) possible, but these are not considered in this thesis. The spacecraft 
kinematics determines the current quaternion state. In this thesis, orbital motion is not 
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considered, so = 0, but it is shown to elarify the relationship between and 
. The rate of the body with respect to the orbit is utilized within the spaeeeraft 
kinematics. This is determined by subtraeting the orbital motion ^(O^ from . 


W = V +V 


W = W-w 


( 1 ) 


Sinee orbital motion is not eonsidered, ^(O^ = . The angular rate of the body 

with respeet to inertial spaee is thus defined in Equation (2). 


w = 


ft). 


( 2 ) 


B, QUATERNION KINEMATICS 

Quaternion kinematies ean be written as a set of four differential equations. The 
first three quaternions ean be written as , while is kept as a separate dynamic 
equation for simplieity, shown in Equation (3) [6, Algorithm 5.76]. 


1 / n r\^ 

q4=--( (O ) Qn, 



(3) 


where 


123 


^1 ^2 ^3 


. The entire quaternion veetor ean be written as 


<1 = 


fl23 




. Another thing to note about quaternions is that the two-norm is 


always equal to one. This allows one to solve for q^, given . The origin for this thesis 
will be defined as^ = [o 0 0 l]- 


11^112==1 

2 2 2 

^4=Vl“^l -^2 -^3 


(4) 
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c. 


MODELING SPACECRAFT KINEMATICS 


Equation (3) is modeled into the Spaeecraft Kinematies Simulink block, shown in 
Figure 3. The direction cosine matrix (DCM) mentioned is only used to determine the 
disturbance torques and is not used here as the disturbances are taken to be zero. 


Spacecraft Kinematics 


dq = 1/2 (q^ x q) 

dq^ = 1/2 ( 0 ( 0 ®)^ q 


Create 

DCM 



D. SPACECRAFT DYNAMICS 


For a rigid body spacecraft with no external torques applied, the total angular 
momentum must be conserved. Since the total angular momentum must be constant, 
then a positive variation in the MED assembly angular momentum will have an 
opposite impact on spacecraft angular momentum and vice versa. All momenta are 
written with respect to the inertial frame. 


hl=hl + h 


where = constant 


( 5 ) 
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Angular momentum is equivalent to the produet of inertia and angular rate 
{h = J(o). Applied speeifieally to the spaeeeraft, represents the total struetural 

inertia of the spaeeeraft. The inertia of the MED assembly is ineluded within the 
spaeeeraft inertia. All configuration changes of the MED assembly with respect to the 
relatively small spacecraft body angular velocity will be considered negligible 
{jmedA ~ O). The angular rate of the spacecraft is expressed in the body frame relative 

to the inertial frame . 

^s/c Jslc^slc s/cbody JmedA ) ^ 


In this thesis, equations often are written in dyadic form to ensure that notation of 
the bases is clearly understood. Some derivations involving multiple frames can be easier 
manipulated (and understood) when utilizing dyadic form. Eollowing this notation and 
assuming a diagonal inertia tensor, Equation (6) can be rewritten as 


Kc = [jsic.kk + jsic,,kk+jy/cJA ) • ((oA+(Oyk+ ) 

= JsicAA [k'k)+JsicAyk [k'k)+JsicAzk [k'k] 

+ Jsic^Azck [k • )+JsicAyk [k • 4 )+JsicAzk [k • k ) 
+ isic^Azck [k • )+isicPyk [k • 4)+ isi^^zk [k • k ) 


(V) 


Since are orthogonal, S, = 1 and h, •bj=0 . This leads to 


Kc = A/c. + j.dc,Ay^2 + j.dc,Az^3 • 


( 8 ) 


Sometimes it will be more convenient to show an equation in matrix form. 
Equation (9), but either form can be used universally. 



J sic^^ 

0 

0 

B 


B 

Js/cA- 

II 

0 

J si^22 

0 




Jslc22^y 


0 

0 

J 5/C33 


^z 


js/c,,^Z 


( 9 ) 
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Since angular momentum must be constant, the derivative with respect to time 
must be zero. 


=0 


( 10 ) 


External torques ean be added as disturbance torques per Equation (11). 




( 11 ) 


Disturbanee examples in EEO are gravity gradient, aerodynamie loads, and solar 
radiation. These and other external torques will be eonsidered negligible in this thesis and 
will be set to zero (T^, = 0). 


Equation (10) reveals that the spaeeeraft and MED assembly angular momentum 
rates must be equivalent in magnitude, opposite in direction. Also, the torque required to 
manipulate the spacecraft corresponds to the spacecraft angular momentum derivative in 
the inertial frame . Thus, the two terms ean be used interchangeably to express 

the torque applied to the spacecraft by the MED assembly. Aecordingly, the torque 
required of the MED assembly is equivalent to the negative torque in the body frame, 
shown in Equation (12). 




= -h 


sic 



( 12 ) 


To obtain the angular momentum rate with respeet to the referenee frame, the 
angular rate between the frames is necessary. Therefore, the angular momentum rate 
(torque) with respect to the inertial frame is equal to the sum of the torque in the body 
frame and the eross produet between angular rate and angular momentum. 




(13) 


The time derivative of Equation (6) is shown in Equation (14). The 
eomponent is only applicable if fuel is being expended during the maneuver or if flexible 
elements are causing change. It is assumed that any ehanges would be very small; 
therefore the spacecraft can be modeled as a rigid body. Moreover, it is assumed that the 
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inertia of the spaceeraft is not changing over the course of a maneuver and the term is 
zero, shown in Equation (15). 


j , j N -B 

h^r ~ J dr ^ J dr ^ 


UB J N -B 

^s/c ~ J s/c ® 


Substituting Equation (6) and Equation (15) into Equation (13) yields 


|1 V J N-B , N^^B . . J N^^B 

hdr ~ J dr ^ “I" ^ ^ J dr ^ 


Equation (17) shows the derivation in matrix form. 


hl = 


Jslcn^x + (i)/c33 “ hlC22 )^z^y 
jslC22^y + (i)/c„ “ 7 )/C 33 
Js/C22^Z + {js/C22 ~ hiCu X 


(14) 

(15) 


(16) 


(17) 


The derivation may also be carried out in dyadic form: 

Kc = hlc„(^xk + jslC22^yk + jslC22(^zk 

(18) 

+ (®A + ® A + O^zk)^ (a,.® A + js/C22^yk + js/C22^zk) 

k/c = jslc,,(oA + jslC22^yk + jslC22(^zk 

+(^33 - )«,«A+(a. - A 33 +(A 3 - A.aA 

= (a Ax+(As “ A 2 ) ® A. )^i ( 19 ) 

+ ( Ac33®,. + ( - Ac33 )®X®Z )4 
+ ( Ac33®Z + [islC22 - isle,, )" A^ ) A 


Correspondingly, by referring back to Equation (12), the torque of the MED 
assembly is related to the spacecraft torque per Equation (20). The MED portion is 
written in the inertial frame and includes the stored momentum in the body frame, to be 
later derived in Chapters IV and VI for reaction wheels and CMGs, respectively. 

-i\N J N \B N B . - T N B /'^fW 

^„.x=~Jsic (O - m Q) (20) 
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Rewriting the angular momentum rate in Equation (20) as its eorresponding 
torque yields 



medA 


-T - (O XJ , CO 


( 21 ) 


E, MODELING SPACECRAFT DYNAMICS 

As a slight variation of Equation (21), Equation (22) is modeled within the 
Spaeeeraft Dynamies Simulink bloek, shown in Eigure 4. 


_B _]V N -B .. j N ..B 

T =-T - (O XJ, (O 

sic medA 




Note that the disturbanee torques = 0 . Thus, this term is not ineluded in (22). 


Spacecraft Dynamics 



Eigure 4. Spaeeeraft Dynamies Model 
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III. EIGENAXIS REORIENTATION 


While traveling at several thousand kilometers per seeond, a spaeeeraft is 
eommanded to move from one quaternion state to another. This reorientation time may 
take a second or two for small maneuvers or tens of seconds to minutes for large 
maneuvers. For LEO Earth-imaging spacecraft, this reorientation must occur swiftly and 
precisely. Reorientation (slew) time is time lost as it is time not conducting mission- 
essential targeting. Since spacecraft can traverse key areas in minutes (or possibly 
seconds), slew time ultimately limits the number of targets per pass. Since these multi¬ 
million or billion dollar spacecraft are only productive when they are not slewing, slew 
time is productivity lost. 

A. SHORTEST ANGLE MANEUVER 

A simple, effective slew maneuver is the eigenaxis reorientation or eigenaxis 
slew. The eigenaxis slew is a rotation along the shortest arc between two attitudes. This 
maneuver is used due to its simplicity, predictability, and reasonably rapid speed, as 
opposed to performing successive rotations about the individual body axes. Given three- 
axis control, there are infinite paths that may be taken to transition between two attitudes. 
The eigenaxis slew utilizes the shortest of these paths. The axis of rotation is known as 
the eigenaxis, whose direction remains constant and orthogonal to the direction of 
angular travel throughout the maneuver. Assuming that the spacecraft can be reoriented 
while maintaining rotation about the eigenaxis, the problem can be simplified into a 
single angle slew, shown in Eigure 5. 


17 




eigenaxis 
(out of page) 


6 rotation angle 


Figure 5. Eigenaxis Slew 


For a reorientation from one angle to another, a direction cosine matrix (DCM) 
(discussed in greater detail in Chapter IV) can be used to show the necessary 
relative reorientation from reference A to reference B. There is at least one eigenvector 
associated with the DCM that produces a unity eigenvalue shown by Sidi [5, Algorithm 


4.1]. 


R^e = le 


B 


(23) 


Equation (23) represents the fact that instead of breaking a DCM into a series of 
three consecutive rotations, the reorientation can be created by using a single rotation 
about a specific axis or eigenvector. Unity eigenvalue occurs when rotating about the 
principle eigenvector or eigenaxis. Because of this, the eigenaxis slew is the shortest 
angle maneuver. 

Assume the spacecraft’s orientation is originally described by quaternions such 
that f with a requirement to slew to some final quaternion 

state = [ qf^ f . Note, q^ is the scalar. The quaternion difference q^ 

can be described by Equation (24) [6, Algorithm 7.69]. 
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(24) 





9^/3 
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-9'/3 
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9^/3 


00, 


Equation (25) shows the relationship between the quaternion difference, the 
eigenaxis, and the slew angle [5, Algorithm 4.11]. 


=esin 




cos 


A 

V 2 y 


(25) 


\ ^ J 


Rearranging Equation (25), the desired slew angle can now be solved for. 

0, = 2cos-‘(^^J 


(26) 


The eigenaxis can now be solved as 

1 

sin(0^/2)^'"'^^ 


(27) 


B. QUATERNION FEEDBACK 

Using an eigenaxis slew, the reorientation problem can be simplified into a single 
angle slew, like the tracing of an arc of a circle on a 2D plane. This section will explain 
how the complicated spacecraft dynamics derived in Chapter II can be modeled as a 
linear double integrator, providing the spacecraft travels precisely along this eigenaxis 
arc. The double integrator model will be obtained from the eigenaxis quaternion feedback 
(EQE) control logic often used in spacecraft ACS [6, Algorithm 7.74]. 

«" = -kJs,c 9,,3 X J.c , (28) 

Recall the nonlinear spacecraft torque from Equation (16). This inertial torque is 
what the EQE controller is attempting to control [u" = . 
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, 11 J N'B , J N^^B 

« =T^s/c = Js/c (0+0) a 


(29) 


Equations (28) and (29) both contain the same gyroseopie eoupling term, a term 
not eondueive to solving a linear trajeetory along an are. The nonlinearity ean be 
eliminated and a new eontrol u’ is determined in Equation (30). This shows that the 
gyroseopie eoupling term was simply added to the EQE to eliminate this nonlinearity in 
the eontrol. 


= (30) 


Sinee = JO)’^ from Equation (15), the eontrol u’ is effeetively a torque 
eontrol in the body frame. Equation (30) ean be rearranged as follows. 


«' +kj,,q,^^^+cj,/co^=0 

J,/d)^ + kJ,^q^^^^+cJ,/a)‘^=0 


(31) 


The eonstant inertia ean then be divided out to assume aeceleration as the eontrol 
variable. 




(32) 


u +kq + c ^ co^ — 0 

-*^123 

N • B , 1 ^ I N ^^B r\ 

CO + k q^ +c CO =0 

-*^123 


(33) 


The problem thus far is measured in three-dimensional space. But the purpose of 
an EQE is to maintain the trajectory along the eigenaxis. Sinee the arc is measured along 
a plane, the quaternion trajeetory ean be measured approximately as a single angle, while 
the body rate ean be measured as a single angular rate, both measured as an error. The 
simplifieation is explained in [6, Example 7.15]. 

d, = d-d, (34) 

u +k{0-0^) + c(0-0,) = 0 

e, + k{e-e,)+c(e-e,) = o 
20 


(35) 



In Equation (33), ^(O^ was not considered an error term because the control 
assumes no shaped angular velocity. Instead, it regulates it to zero. Therefore, the rate 
error is equal to body rate 0^ = 0= . The fundamental formula for a proportional 

plus derivative controller (PD) is shown as 

0£ + ^0£ + c0£=O (36) 


Equation (36) is now written in terms of a double integrator, where the control is 
equal to the second derivative of the signal u = 6^ = a . 


u = a 

= -kd^-cd^ 


(37) 


Assume a desired angle and rate are given. The controller can now regulate the 
acceleration signal to drive the error to zero. As stated previously, the desired rate is 
zero, therefore it is eliminated from the equation. This means that Equation (28) was 
treated as a rest-to-rest maneuver. The double integrator in Equation (38) is implemented 
in Figure 6. 

u = a = -k[0-0^)-cco (38) 


Desired Input Controller Spacecraft 

Input Error ^^ ^_|_ 




Figure 6. PD Control of a Double Integrator Model 

Therefore, the EQF controller is essentially just a PD control system in disguise. 
The gyroscopic coupling term can be added back in to the PD controller, just to be 
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subtracted back out by the spacecraft physics. Therefore, the spacecraft can be slewed as 
a double integrator, as shown in the simplistic spacecraft simulation in Figure 6. 


C. GAIN DESIGN FOR QUATERNION ERROR FEEDBACK 

If the goal is to perform a rest-to-rest reorientation, may be fixed and act as a 
step input, which moves from zero to the desired angle. As with any step input, the error 
is very large at first and decreases significantly over time. Proportional and Differential 
Gains {k,c) must therefore be chosen carefully. 


Underdamped systems are the quickest to respond, but must be carefully designed 
so the torque (acceleration) and overshoot do not become too large. The percent 
overshoot (Mp) is a measure of the overshoot compared to the starting/ending condition 
difference. For the following example, 10% overshoot is used. Settling time is the time 
it takes the system to fall and remain within a certain band, commonly determined by 
imaging and targeting requirements (e.g., 2% of the final or steady-state value). An 
example of gain selection is shown in (39) such that the percent overshoot and settling 
time is the selection criteria for damping ratio and natural frequency [22, Algorithms 
4.39, 4.42]. 


c 


-ln(Afp/100) 
^;r" + ln(Mp/100) 


CO 


n 



(39) 


The previously mentioned k and c are the feedback gains for the second order, 
homogeneous system shown in Equation (40). The same canonical expression is shown 
in Equation (41) in terms of the damping ratio and natural frequency. 


e+ce+ke=0 

(40) 

e + 2^coj + co^^e = 0 

(41) 
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Now, it can be seen that 


k = a,^ 


(42) 


Using the feedback control system mentioned in Section C, the ideal gains for a 
quick 10° slew will be much different than those for a 60° slew due to the fact that the 
control system operates proportional to the angle error, given the same hardware. As an 
example. Figure 7 shows the time history of the system variables for the underdamped 
case. Both 10° (solid lines) and 20° (dashed lines) reorientations occurred using the same 
set of gains. Notice that they both settle at the same time. The 20° slew, however, shows 
significantly higher acceleration and rate commands. The phase space indicates that both 
slews follow essentially the same trajectory, and the magnitudes are scaled since the 
system is linear. 



Phase Space 


10 

time (sec) 



Figure 7. Underdamped Reorientation (solid 10°, dashed 20° slews) 


It would be expected that a 10° slew could occur about twice as fast as a 20° slew 
for a given value of peak acceleration, but naturally, the linear control system prevents 
this. As a further example, the same 10° slew is compared with a 60° slew to show 
further exaggeration of this restriction. Figure 8 demonstrates that the required 
acceleration and rates increase significantly with the slew size. 
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Control and State Variable Trajectories 


Phase Space 




0 

Figure 8. 


5 10 15 20 0 20 40 60 80 

time (sec) 9(t) 

Underdamped Reorientation (solid 10°, dashed 60° slews) 


At some point, the aetuators will not be able to produee the required torque 
(aeceleration) and/or angular momentum (rate). Suppose the 10/20° slews from Figure 7 
are now subjeet to torque and momentum saturation limits that are equivalent to 5°/s^ and 
4°/s, respectively (values representing a very agile spacecraft). This means the MED 
array cannot produce more than these limitations around the desired eigenaxis. Figure 9 
demonstrates the saturations by flat regions, where maximum torque and acceleration are 
held for extended periods of time. The phase space indicates both trajectories were 
saturated, but the 20° slew was saturated twice. Within the same 20 seconds, the 20° slew 
was not complete due to the torque and momentum limits. 
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Phase Space 


Control and State Variable Trajectories 




0 5 10 15 20 0 10 20 30 

time (sec) 9(t) 

Figure 9. Underdamped Reorientation (solid 10°, dashed 20° slews) a, ro Limited 


Notice that the spacecraft began at maximum torque and held it until the 
momentum limit was reached. At that point, maximum momentum was held until 
crossing the target angle, after which the spacecraft decelerated at maximum torque. 

Engineers typically want to optimize trajectories, minimizing overshoot and 
settling time to reduce wasted transient distance and time. To perform a slew quickly 
with no overshoot, the critically damped case may be tried (Figure 10). 


Control and State Variable Trajectories 
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Figure 10. Critically Damped Reorientation (solid 10°, dashed 20° slews) 
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Notice that in Figure 10, both slews do not require anywhere near maximum torque 
or momentum. In this simulation, using the eritieally damped gains seems like a good 
approaeh. In a real spaeeeraft, the mass (inertia) will deerease over time as fuel is 
expended. Thus, spaee system engineers must plan in a buffer factor in order to account 
for the ehange in mass. As a “rule of thumb,” an engineer may select an initial damping 
ratio around ^ = 0.7. As mass decreases, the damping ratio inereases since the two 

values are inversely proportional 1 / m). The damping ratio will eventually approach 

critical damping and move into overdamped towards the end of the life of the spaeeeraft. 
Thus, one eould argue that the eritieal damping or near-eritieal damping seenarios would 
not be useful in spaeeeraft design. 

A simple approaeh to eontrol system design is to devise a set of gains that work 
aeross a large range of slew angles. To satisfy all angles, the eontrol engineer must design 
the gains for the largest slew angle the spaeeeraft is regularly expeeted to perform. If this 
is done using a linear model, the system may perform very poorly for small angle 
reorientations. Therefore, the engineer ean be foreed to overdesign the hardware to satisfy 
the torque and momentum requirements for smaller angle maneuvers. 

D, TRAJECTORY SHAPING 

Gain design is very delieate beeause the engineer must tune the gains speeial for 
the spaeeeraft, regardless of slew size. Instead, input shaping allows a fixed set of gains 
that must eorrect for only a diserete minimal error to be utilized for multiple slew sizes. 
Optimal trajectories ean be pursued without major hardware upgrades. 

Input shaping is an alternative to using a step eommand as the input to the attitude 
eontrol system. Input shaping ean be introdueed in the form of a set of shaped torque 
eommands, represented as an aeeeleration input in Figure 13 [17]-[18]. This is ideal, 
beeause it completely avoids the need for a eontroller. Torque/aeeeleration input shaping 
is, however, unrealistic because the system would be operated in the open loop. This 
means that it does not have a feedbaek meehanism to aeeount for uneertainties in the 
spaeeeraft inertia and external disturbanees on the spaeeeraft, whieh is undesirable. 
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Figure 11. Torque Input Shaping Model 


A better approach is to shape the trajectory of the input to a closed-loop system; 
the input angle itself Rather than simply inserting a desired angle 6^ and letting the 
onboard controller calculate how to respond. The desired spacecraft response can be 
decided ahead of time and then implemented as a series of small steps. For example, a 
10° slew can be divided into 100 or 1,000 incremental steps. Now, the gains can be 
optimized for the step size (fixed or extremely small) rather than the entire maneuver size 
(variable). 
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Figure 12. Angle Input Shaping Model 


As long as the spacecraft model is reasonably accurate, the path can be 
predetermined accurately. The feedback loop then adjusts the path to accommodate 
spacecraft uncertainties and external disturbances. Similar methods have been used in at 
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least a few cases in the industry [12]-[14]. The input shaped trajectory 0 ^ is fed to the 
system much like that in Figure 12. 


The ideal input is the trapezoidal input (referring to shape of W ) shown in Figure 
13. This input accounts for both the torque and momentum limits a priori to control the 
error 6 ^ so that the control system does not saturate. 


/' 

time angle 6, 

Figure 13. Rate Limited Eigenaxis Slew 

Due to the shape of the angular rate trajectory, the torque/momentum limited 
input is also referred to as a bang-coast-bang maneuver. This is exactly the command¬ 
generated trajectory used for the Clementine spacecraft that was used for lunar mapping 
in 1994 [12, Figure 2]. The eigenaxis and trajectory angle were first calculated by the 
onboard computer to transit from the current quaternion to the desired quaternion. In 
Chapter V, it will be confirmed that this trapezoidal input is in fact the time-optimal 
solution for eigenaxis slews. 
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IV. REACTION WHEELS 


Reaction wheels have their physical limitations. They have maximum angular 
rates and accelerations that must be modeled. This will ensure that even if a large torque 
is demanded, only the maximum MED assembly torque will be applied. In order to apply 
the saturations, a model of the reaction wheels must be created. 

A basic model of a reaction wheel is shown in Figure 14. The reaction wheel is 
shown in its own local frame, referenced as the wheel frame with its spin axis aligned 
with the m >3 axis. A motor, attached to the body frame, drives rotation of the wheel. The 
principle of operation is for the motor to apply a torque to the wheel. This causes the 
wheel to rotate quickly. The resultant equal and opposite torque is applied back through 
the motor to the body frame. This causes a rotation of the much larger spacecraft body. 
The momentum and torque vectors align with the axis. 





Figure 14. Reaction Wheel Schematic 
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A. 


ORIENTATION 


The angular momentum of a single reaetion wheel (RW) within its (ith) wheel 

frame W with respeet to body ean be shown in the spaeeeraft body frame B through a 
simple rotation. The fundamental formula for this ehange of basis is given as 

h‘ (43) 


In dyadie form, (43) ean be looked at as the projeetion of the wheel frame onto 
the body frame (eomponent-wise dot produet). 





- 











h;^‘ 

(44) 


4 * ^(2 

^ 3 * ^(3 




This transformation ean be eompleted as a matter of three fundamental rotations, 
and represented as a direetion eosine matrix (DCM). A DCM is an invertible square 
matrix where the transpose is the inverse sueh that 

= or J=/ (45) 

This property allows a mueh more straightforward approaeh of first plaeing the 
RW in the B frame, eonverting to the frame via , and then using the transpose 
to eonvert baek. In the example shown in Figure 15a, the MED is plaeed direetly on 
the B xyz origin, rotating about the axis per Figure 14. This produees an angular 
momentum direetly in the ^3 direetion for a reaetion wheel. 
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Figure 15. Reaction Wheel Transformation 


First, a fundamental 3 rotation about the shared axis must be completed to 

orient the w, and axes for subsequent rotations (Figure 15b). This will allow for 
multiple reaction wheels placed along the same skew angle (j8), forming a pyramid. This 
first rotation can be represented in (46). 




cos a,, sina^ 0 

-sino!; cos a, 0 
0 0 1 


(46) 


Assuming a = constant for the initial 3 rotation of the representative wheel, a 
2 rotation is now completed about the axis (Figure 15c). This rotation is used for both 
reaction wheels and CMGs, and is referred to as the skew angle. The skew angle is 
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normally assumed the same aeross all deviees, so the subscript “z” will not be used to 
distinguish which MED it refers to. This rotation can be represented in Equation (47). 




^ cos/3 0 -sin/3 ^ 
0 1 0 

sin 0 cos /3 


(47) 


The previous series of rotations can now be combined into a single transformation 
converting a vector from the body frame to the wheel frame. 
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_1 


ca. 

sa, 
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-sa, 

ca, 

0 

_1 

0 

cP _ 

2 

0 

0 

1 


caffi 

-sa. 

ca^sP 


saf^ -sfi 
ca. 0 
sa^sP cj8 


(48) 


It follows from Equation (45) that '^R^‘ is now the transpose, which is equivalent 
to the reverse sequence via the negative angle of each rotation. Equation (49) is the 
physical representation of the dot product matrix shown in Equation (44). 




B 




-sa. casP 
sa^P ca. sasf5 


-sfi 0 cj8 


(49) 


Note also that 


- 



B 








ca,cj5 

-sa, 

casf5 


k*^,2 

k*w,^ 
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sa,cP 

ca. 

sasP 

1 

• 


k*^i3 


-sf5 

0 

cP 


( 50 ) 
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As an example, plaee the initial MED at angle a = 0°. The resultant 

transformation from the W frame to the B frame is shown in Equation (51) for a 
reaction wheel. 


B nW 


B 


R,{-0°)R,(-I5) = 


cj8 0 
0 1 0 
0 cp 


(51) 


Recall that the RWs in Equation (51) were stationed at a = 0°. A classic 
configuration is using four MEDs in a square pyramid configuration in which each of the 
MEDs are placed 90° apart [15]. Arranged at the proper skew angle, this tricetrix 
arrangement maximizes torques in all directions, minimizing torque cancellations. This 
also provides redundancy in case of a single wheel failure. The arrangement used will be: 
ctj = 0 °,a 2 = 90°,a3 = 180°, 0^4 = 270°. Figure 16 demonstrates this common placement. 





RW4 





Figure 16. Body Frame to MED Rotations 
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Based on the wheel axis description mentioned above, the angular momentum in 
the wheel frames can be written as 


h 


Wi 

nVj 


0 

0 

1 


iH'; 


(52) 


Transforming angular momentum from the W frame to the B frame using 
Equation (43) and « = 0° yields (53). 


or in general form 



sf5 


\B 



hi = 


4*^,3 
K • w. 


(53) 


(54) 


The total angular momentum of the system is simply the sum of each of the 
MEDs in the body frame. 


h 


B 

meclA 



(55) 


Eor the reaction wheel case in Equation (54), the individual rotation matrix and 
angular momentum components can be combined into the product of two matrices, 
equivalent to Equation (55). The first matrix represents the orientation of each wheel with 
respect to the body frame, separated as individual columns. Eor n MEDs, this 3xn 
orientation matrix is multiplied by an nxl vector containing the momentum in the wheel 
frame as shown in Equation (56). 
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r- -| 




^1*^„3 


h, 

(=1 

^2*^13 

4 • ^2, 



hi 


1 

• 

4*^23 

• 

1_ 


_ K _ 


( 56 ) 


The 3xn matrix that relates n MEDs to the principle body axes will be labeled this 
point forward as the orientation matrix . Since it contains only elements from the 3^^^ 
column of the original rotation matrix, it will be further described with a subscript 3: 

Br^W 

“^3 • 




3 


B 


^l*^l3 

• ^23 

brw, 



yV 


K • ^23 




y,. 

^•^l3 

^•^23 



Now, Equation (56) can be rewritten as 




where = 


K 

hi 


(57) 


(58) 


(59) 


For the square pyramid mentioned, the orientation matrix is shown in Equation 
(60) for reaction wheels. 
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wW 



0 

cp 


0 

sP 

cl5 


-sp 0 

0 -sp 


(60) 


35 



Any RW orientation configuration Z can be modeled by specifying 
the appropriate a and P rotations. Note that Equation (60) assumes 

a = [ 0 90 180 270 J. 

Other than the case where n=3, Z will always be a 3xn matrix where n>3. This is 
caused because the system is over-actuated, i.e., more wheels exist than are necessary to 
control the spacecraft attitude. A common approach for performing control allocation in 
the over-actuated case is to use the pseudo-inverse to map three orthogonal commands in 
the body frame to n commands for the individual wheels. 

Z^ = Z^(ZZ^y (61) 


B. DIFFERENTIAL EQUATIONS OF MOTION 

Like the spacecraft, the momentum of the reaction wheel can be defined as the 
product of wheel inertia and angular velocity of the wheel, all with respect to the inertial 
frame. 


ZI — T 

^Wi ~ Jwi 


(62) 


Previously, the wheel inertia was combined into the spacecraft inertia in Equation 

(6), where h^ic-{^Jskbody'^JmedA)^^^- Therefore, the wheel inertia rotating with the 

spacecraft body must be pulled out, leaving only the configuration of the wheel with 
respect to the body. 




(63) 
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( 64 ) 


Shown in dyadic form: 
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( 65 ) 


hj' = )• ) 


Change of basis to the body frame: 


hj = j„i,PiW,^ • (bA + kk + ^3^3) 

= Li,P-i[^i, 'k]k + Li,P-i[^i, 'k)k + 'k)k 


( 66 ) 


The dot produet of two veetor eomponents ean be flipped interehangeably (i.e., 
= w^^ •by). This is important to note beeause the transformation matrix ean now be 
written in the more eonventional notation shown in Equation (44). The resultant rotation 
matrix is shown in Equation (67). 
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b. • w. 
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1 I3 


0 

hJ = 

k*^ly 

k*^i. 

k*^i. 
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1 Q, 



^ 3*^,2 

^ 3*^,3 


tVt33 1 




(67) 


Equation (67) ean be simplified into a version eontaining only the rotation matrix 
portion multiplied by its momentum, expressed in the body frame. 


hj = 


k'^i, 

k'^i. 


/ Q 

J H’t33 I 


( 68 ) 


Erom Equation (56), the total angular momentum for n wheels in the body frame 
can be rewritten as 
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K * 1^23 

4*^23 




Jw ^1 



./ Wj 33 1 

1 /«3 

yv 
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2^233^^ 

^2*W„3 





/ Q 

J w„ 33 n 


( 69 ) 


Equation (69) is equivalent to Equation (58), where hi = . Within the 

model, the inertia and rates are expanded out in the following manner to keep traek of 
them separately. 
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(70) 


The torque with respect to the wheel frame is the time derivative of Equation (70). 


ly. 


= ■ 


dt 


{h:)=j..n’' 


( 71 ) 


The torque with respect to the body will be defined as the change in angular 
momentum with between the wheel and body frames. 


T® = where hi = h, 




(72) 


Expanding out these terms following the chain rule yields Equation (73). Since 
this is a reaction wheel, the orientation is fixed and therefore = 0. The inertia is 
assumed to not change with time, J^„ = 0 . 



B^W T 




(73) 

(74) 


In dyadic form, a single reaction wheel torque can be found using the derivative 
of Equation (66). 
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(^,3 * )^1 + (^<3 * ^2 )k + Li,Pi (^,3 * ^ 3)4 

+ J™33^<(^<3 *^1)^1 +^33^<(^-3 *4)^2 +^,-33^<(^.'3 ’ ^3 ) ^ (75) 

+ j^i,Pi (^,3 * ^1 )^i + iwi,Pi (^,3 * ^2 ) 4 + j«i,Pi (^,3 *^3)4 

The inertia of the wheel is not ehanging with respeet to the wheel frame. Also, the 
orientation of the 3 axis of wheel frame is not moving with respeet to the body frame. 


L, =0 andw. =0 


Equation (75) now beeomes 


-[hj ) = (w. •b,)b,+ (w,^ •b^)b^ + (w. • 4)4 


Jw, 33 ^^,( 4 *^, 3 ) 

'-(hj)- pA(^2-p) 


brw. 


j^iPi 




By Euler’s transport theorem adopted from [23, Algorithm 8.18a], the inertial 
derivative of Equation (74) is 


’d(h:)_^d(K) 


+ Wx/lf 




Equation (79) ean be rewritten by substituting the body torque from Equation (74) 
and body angular momentum terms as 


_V BryW J rt 

+ CO X J^Q 


By substituting Equation (80) into Equation (12), knowledge that torque is in the 
wheel frame as defined in Equation (71), and a rearrangement of terms, 

the result is an equation in terms of the torque of the reaetion wheel. This equation now 
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involves the pseudo-inverse to translate torque commands given in the body frame to n 
wheels, shown in Equation (82). 

_N BryW -W . N ^ BryW J /QIN 

V BryW~\'^ I _V N ..B BryW J /-J IT \ ^o^'v 

= L ■^3 J [-^re, - (O ^ J^Q j (82) 

Further, Equation (82) can be reduced to Equation (83), solving for the reaction 
wheel dynamics in terms of the spacecraft. 

^ =J« L ■^3 J [-^req - (O ^ j (83) 

C. MODELING REACTION WHEELS 

In order to apply the physical limits of the reaction wheel actuator, the torque 
required of the assembly must be broken down, torque and wheel speed limit saturations 
applied, and then the torque of the reaction wheel assembly rebuilt in Figure 17. Equation 
(82) is used to break down the torque required and (80) is used to rebuild it. 


Z (3xn) 




Z rw 


-KEO 


Reaction Wheel Dynamics 

(saturation occurs in red) 


T^ = Z*(.J^ .N^BxZhW) 
rw rw req rw 


7^ =Z + 

rw rw rw rw 


N_B 


(w_BN)-<«^ 




Cross 


Matrix 

Product 




Cross 
Product h 


(Tq_req) -lg%(^ 


.w 


hw 


rw“ 




.w 


rw 


1/s 


ED^ pinv(Z) 


Matrix 

Multiply 


rW 

rwi 




w 


rw“ 


(^)-( ^Tq_rwA) 


E> Matrix 


Multiply 


Figure 17. Reaction Wheel Model 


This model is contingent that momentum in the wheel frame is derived directly 
from wheel torque from Equation (71) such that 
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h'" = . 

w j*'w 


( 84 ) 


The wheel torque limitations oeeur in the individual wheel frames. The rate 
limitations are modeled as wheel momentum saturations. Note, this model does not shut 
off the wheel torque when ||/j,||=• It is not neeessary for the eomparisons in this 
thesis beeause only torque and momentum limited signals are used. 


D. REACTION WHEEL MOMENTUM SPACE 

By measuring the maximum momentum in any given direetion over a In 
steradian for a partieular reaetion wheel eonfiguration, a boundary ean be defined that 
deseribes the maximum momentum eapaeity of the array. This is the physieal limitation 
of the system defined by the veetor sum of the maximum eapaeities of the individual 
momentum wheels. If this boundary is ereated for every possible direetion around a unit 
eirele in the body frame, the surfaee of boundaries is typieally known as the momentum 
envelope [15, p. 1607]. To determine the best possible meehanieal eonfiguration of 
systems, the volume must be maximized, ideally forming the most spherieal shape 
possible. Sinee eaeh of four wheels is maximized individually, the resultant shape has 
radii larger than . If all four wheels were aligned, a maximum radius of 4 eould be 
aehieved, leaving zero momentum in other direetions. The maximum momentum spaee 
for an evenly-spaeed (Act = 90°) four wheel system is aehieved when 


P = cos ' 


(—] 

vV3. 


raJ = 54.74° . 


(85) 


The maximum angle in Equation (85) is shown in [15, p. 1610] as the conjugate 
35.26°. The resulting envelope is shown in Figure 18. The figure was created using a 
mesh of triangles between all surface points. 

Assuming all four wheels have equal momentum storage. Figure 17 is calculated 
as a function of the maximum momentum available to any given wheel . The 
envelope scales symmetrically with the size of the wheel ). Therefore, larger wheels 
would push the surface outward, allowing more momentum capacity, but the shape would 
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be retained. This is why the dodecahedron is generic for any four-wheel system, but 
unique to only four-wheel systems. A five-wheel system would have a much different 
shape with a maximum at a different angle {3 and a larger volume (see [15]). 
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Figure 18. Momentum Space for 4 Wheel Reaction Wheel (a= 0,90,180,270°) 


The momentum envelope shown represents the angular momentum available 
about any given rotation axis. During times of maximum momentum use, the momentum 
vector would touch the outside surface of the envelope until the system begins to 
decelerate and momentum decreases back to the center of the shape where momentum is 
zero. The vertices of the momentum envelope indicate optimal axes that allow locally- 
rapid rotation rates. Although fourteen vertices exist, only six are global maxima (marked 
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with blue). The centers of each facet indicate the minima for each respective plane. 
These minima are also the global minima, each marked in magenta. 

The four-wheel polyhedron is again shown in three orthogonal views in Figure 19 
to clarify its shape. The blue maximum vertices (global maxima) occur at magnitudes 

equal to = 2309h^^^. The black minor vertices (local maxima) occur at 

^max “ ■ Thc minimum /i® is actually zero, but represents the minimum across 

the momentum envelope. The minimum radius occurs at the middle of the flat facets 
where = 1 , interestingly equal for every axis. 



Max Vertices - 2.309 /j^ 
Minor Vertices |//| = 2.000/?^^^ 
Minimum |/j| = 1.633/?..,,. 


'^2 

Figure 19. Momentum Space of a 4 Reaction Wheel System 
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Table 1 lists all of the possible eonfigurations of momentum vertiees and faeet 
minimums. Onee normalized, these represent the respeetive axes with the momentum 
magnitudes listed above. 


Max Vertices Minor Vertices Minimum Facets 

[ 1 1 0 ] [ sin(P) 0 cos(P) ] [ 1 0 0 ] 

[ 1-1 0 ] [ sin(P) 0 -cos(|3) ] [ -1 0 0 ] 

[ -1 1 0 ] [ -sin(|3) 0 cos(P) ] [ 0 1 0 ] 

[ -1 -1 0 ] [ -sin(P) 0 -cos(P) ] [ 0 -1 0 ] 

[ 0 0 1 ] [ 0 sin(P) cos(P) ] [ sin(P) sin(P) 2cos(P) ] 

[ 0 0 -1 ] [ 0 sin(P) -cos(|3) ] [ sin(P) sin(P) -2cos(P) ] 

[ 0 -sin(|3) cos(P) ] [ sin(P) -sin(|3) 2cos(|3) ] 

[ 0 -sin(P) -cos(P) ] [ sin(P) -sin(P) -2cos(P) ] 

[ -sin(P) sin(P) 2cos(P) ] 

[ -sin(|3) sin(P) -2cos(P) ] 

[ -sin(|3) -sin(|3) 2cos(|3) ] 

[ -sin(P) -sin(P) -2cos(P) ] 

Table 1. Reaction Wheel Significant Axes 

From Figure 17, the wheel torque vectors directly align with the wheel 
momentum vectors. This causes a torque envelope equal in shape to that in Figures 18 
and 19. Instead of using as the scaling factor, the dodecahedron would be instead 
scaled by the maximum individual wheel torque . All maximum and minimum 
vertices are co-aligned with the momentum vertices, just with different units. 

E. LARGEST INSCRIBED SPHERE 

As mentioned in Chapter I, sometimes it is necessary to create artificial 
limitations on a control system. Recall from Figure 18 that the twelve facet minima occur 
at = 1.633/2^^ . Although they are considered minima, they still lie on the momentum 

surface. For example, on axis[ 1 0 0 ], momentum up to 1.633/2^^^ (torque up to 
1.633 Tj^^^) is available. This is the physical system limit that cannot be exceeded for an 
eigenaxis rotation about that specific axis. Since this is the global envelope minimum, a 
sphere can be created with radius that includes all twelve points. Such a sphere 
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is the largest possible sphere within the momentum space, shown in Figure 20. 
Restricting commanded momentum to lie within the inscribed spherical volume ensures 
that 1.633hj^^^ is available about every possible axis. 



Figure 20. Reaction Wheel Momentum Space with Largest Inscribed Sphere 


Of course, 40% more momentum is available about the maximum axes shown 
in Table 1, but in a conventional attitude control system design, this additional 
torque/momentum capacity is sacrificed. This sacrifice is made to ensure a designed slew 
rate is available for any given slew. The inertia of a spacecraft is rarely symmetric and 
equal about every major axis. The resultant is an inertia ellipsoid, further reducing the 
momentum space envelope. This occurs due to the need to ensure the momentum does 
not change with a non-symmetric inertia matrix. The maximum spacecraft inertia value is 
also used about every axis, further limiting the system to enforce consistent slew rates. 
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F. 


REACTION WHEEL PSEUDO-INVERSE STEERING LOGIC 


The physical limitation of a reaction wheel system was previously discussed. 
Another limitation associated with reaction wheels is the mathematical limitation of the 
Moore-Penrose pseudo-inverse, Equation (61). When reorienting coordinates from the n 
wheels in the wheel frame W to three orthogonal coordinates of the body frame B , there 
is a unique solution. 




( 86 ) 


The pseudo-inverse offers a single solution for translating the actuation from the 
body coordinates (3x1 vector) uniformly across n wheels {nxl vector). The pseudo¬ 
inverse is desirable because it is the minimum least squares (L2) solution, minimizing the 
sum of the squares of all of the wheels’ individual momenta/torque [15, p. 1609]. 


hW^^BzWjh’^ 

^W^^BzWj^B 


(87) 


Since the system is underdetermined, there are infinite possible nxl wheel 
configurations that will satisfy the body coordinate solution, but the pseudo-inverse will 
give only one such configuration. 


Pseudo-Inverse Control Allocation 


Using the reaction wheel orientation matrix for a 4 wheel system 


(60) and (3 = cos ‘ 


^ 1 ^ 


vV3y 

inverse used in this example is 


rad = 54.74°, the resultant orientation matrix 
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and pseudo- 


Br,W 
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- V2 0 

0 -^/2 

1 1 


( 88 ) 
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yfl 0 1 

2^ 0 V2 1 

4 -V2 0 1 

0 -V2 1 


( 89 ) 


To illustrate the implieations of , consider the following possible wheel 

momentum configuration where represents a scalar of the maximum possible 

momentum for a given wheel. 

1 

^possible ~ ^ ^max (90) 

1 

Using Equation (90) leads to the following possible momentum in the body 

frame. 

1.633 
0 

1.155 

Assume this possible h'^ was commanded rather than derived. Now, the pseudo¬ 
inverse would be required to derive the individual wheel momentum commands, . 

Applying the least-squares solution [ leads to a different solution than 

1^ 1 1 -1 1 ]^)- The pseudo-inverse solution is actually physically unachievable 
because the largest momentum command, 1.5 , is higher than the maximum for a 

given wheel: 
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pseudo-unacheivable 
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possible 



1.633 
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1.155 


1.5 

0.5 

-0.5 

0.5 




h 


max 


(92) 


In order to implement the least squares eontrol alloeation, the proportionality of 
the solution must remain intaet [24, pp. 880-883]. The wheel momentum vector in 
Equation (92) is therefore normalized by the highest element (greatest in magnitude). 
This allows appropriate control allocation to rotate the spacecraft predictably about the 
desired axis. 
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(93) 


The normalized solution Equation (93) is far different than the wheel momentum 
configuration at the beginning of this example, Equation (90). Translating back to the 
body frame reveals that the pseudo-inverse limited solution in Equation (94) can only 
achieve 67% of the possible momentum shown in Equation (91). 
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(94) 
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The previous example ean be redone using 



( 95 ) 


with a as anything from 0.333 to 1.0. The system returns the same limit at 

r 1.089 T 

Keudo-X^ = 0 ( 96 ) 

0.770 

whieh is only 67% of the possible body momentum in this example. 

r 

1.633 

^possible ~ 0 ^max (97) 

1.155 

The same performanee limitation results in a loss of torque per Equation (87). 


2, Performance Loss 

The momentum veetor in the previous example eorresponds to a minor vertex (see 
Figure 19). 

g r 

1.633 c[5 

Kossihie^ 0 /z_= 0 2.0/z_ (98) 

1.155 J sp 

The pseudo-inverse limited veetor is 67% of giving a magnitude of 

|A| = 1-33A™. 

g r “I® 

1.089 cp 

Keudo-iun^ 0 0 l-33/z„^, =0.67/1 V,,, (99) 

0.770 J sp 
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Graphically, would be well within the reaction wheel momentum 

envelope. Shown in Figure 21, the red shape represents the full reaction wheel 
momentum space. It is the physical limitation of the reaction wheel system utilizing the 
maximum momentum possible to rotate about any given axis. The black dot represents 
the original point (minor vertex) shown in Equation (98). The red dot represents the 
pseudo-inverse limited system shown in Equation (99). The blue/black lines represent all 
possible magnitudes that share the same eigenaxis as the minor vertex. The blue line 
specifically represents the set of all achievable configurations by the pseudo-inverse. The 
black line represents the set of all unachievable, yet still possible, configurations. 






Unachievable 


Achievable 


Minor Vertex 
Possible Magnitudes 
Pseudo-Inverse Limited Max 

Possible Pseudo-Inverse 
Limited Magnitudes 


^1 


Eigure 21. Reaction Wheel Momentum/Torque Space Showing Pseudo-Inverse Eimit 


The same is true for all other minor vertices, as they are unachievable using the 
pseudo-inverse. On the other hand, the minimum points and the maximum vertices are all 
achievable, shown in Eigures 22 to 23. 
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Minor Vertex 

Possible Magnitudes 

Pseudo-Limited Max 

Possible Pseudo-Limited 
Magnitudes 

Pseudo-Limited Space 


-1 



1 


Figure 22. Reaction Wheel Pseudo-Inverse Limited Space 


The representation shown in Figure 22 makes it clear that the pseudo-inverse is a 
mathematical limitation on the system performance and not a physical system limitation. 
Indeed, this is true because the pseudo-inverse minimizes the sum of the squares of the 
moment/torque in a given direction. The least squares solution discards some of the 
available moment capacity in favor of an L2 control allocation. In reality, the spacecraft 
has the ability to rotate faster (higher momentum) about many axes, but the pseudo¬ 
inverse does not necessarily allow this. Going from wheel space to body space, there is 
only one mathematical answer. Going from body space to wheel space, there are an 
infinite number of potential arrangements. The pseudo-inverse only provides the least- 
squares solution, a mathematically limited solution. 

Figure 23 shows the pseudo-inverse limited momentum space from Figure 22 
overlaid into the full reaction wheel space from Figure 19. This makes it clear how much 
volume is missing across the different orthogonal views. An interesting note is that the 
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pseudo-inverse solution intersects all six of the maximum vertices and all twelve of the 
minimum points, but none of the eight minor vertices. 



Figure 23. Reaction Wheel Pseudo-Inverse Limited Space - Orthogonal 


G. PSEUDO-INVERSE LARGEST INSCRIBED SPHERE 

As previously shown with the full reaction wheel space, the full momentum 
envelope can be further limited to an enclosed sphere (Figure 19). Similarly, the pseudo¬ 
inverse limited space can be confined to a sphere to ensure a constant value of achievable 
momentum in all directions. The minimum pseudo-inverse envelope surfaces occur at 
1.333/Zjj^^. Therefore, the pseudo-inverse space can be confined to a sphere with a radius 
of (Figure 24). This will ensure that this momentum is available in all 

directions. 
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Figure 24. Pseudo-Inverse Momentum Spaee with Largest Sphere 

Figure 25 shows a eomparison of the reaction wheel momentum space under the 
various control limitations. Here it is obvious the huge loss of momentum space that 
occurs under the typical spherical simplification utilized in the design of a spacecraft 
attitude control system. Later chapters will be used to discuss methods, particularly the 
use of optimal control, to improve the usage of the available system capacity. 


53 



















hi 

Figure 25. Momentum Space Comparison 
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V. OPTIMAL REACTION WHEEL ATTITUDE CONTROL 


This chapter will be used to derive the time-optimal slew for a reaction wheel 
system. First, the time-optimal eigenaxis slew will be generated. This standard heritage 
maneuver is best only if the engineer wants the spaeeeraft to stay on the eigenaxis to 
minimize the overall rotation angle. The spaeeeraft often has more momentum and torque 
available to exploit, but it is necessary to deviate from an eigenaxis path in order to 
aecomplish this. The remaining seetions of the ehapter will expand the linear derivation 
from Seetion A into the nonlinear reaction wheel and spaeeeraft dynamics developed in 
Chapters II and IV. DIDO, an optimal control solver [21], will be used to obtain the 
optimal control solution. An attitude control architecture will also be presented in order to 
implement the optimal solution. 

A. TIME-OPTIMAL EIGENAXIS REORIENTATION 

The time-optimal, rest-to-rest eigenaxis slew solution is fairly straightforward. 
Not all of the following steps are necessary to solve this problem, but all steps will be 
covered for completeness. The optimal control methodology is necessary to understand 
prior to Seetion B, where a mueh more eomplieated derivation is presented. 


1, Problem Formulation 


In general, figures of merit are used to determine how well the slew meets the 
requirements important to the engineers and/or customer. These are broken up into eosts 
(stuff to minimize) and eonstraints (boundary conditions that must not be exceeded). The 
maneuver is defined as rest-to-rest with the initial and final angles predetermined. 

Initial time is zero and the final time is free (term to be optimized), if it is desired 
to design a rapid slew maneuver. 


t, = f=0 


tf = V 


( 100 ) 
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As shown in Chapter III, the spaeecraft model can be simplified as a double 
integrator. In this case, the state variables and control are 


1 

_ 1 


1 

1 _ 

1 

1 _ 


1 

1 _ 


x(t)- 

U = [m( 0] = [ct(t)] 

The state space representation of Equation (101) is shown as 


( 101 ) 

( 102 ) 


x = f{x(t)Mt)) 


1 

_1 


1 

_1 


1 

_1 

1 

1_ 


1 

1_ 


1 

_1 


(103) 


First and foremost, the solution must meet the boundary conditions. The purpose 
of slewing is to accurately point at a new target (end condition). If the spacecraft slews 
but does not get the desired end condition, then it may not be possible to perform the 
mission requirements. The final angle will be defined as the difference from start to 
finish. Therefore, the initial angle will be set as zero. Initial and final rates are zero 
because the maneuver is considered rest-to-rest. The boundary conditions are defined as 
follows. 


^(C) = ^o =0 

O)(to) = 0)^=0 
O(tf) = Of=0j 
= =0 


These boundary conditions can be rewritten as 


x(to) = A:° = 


x(t^ )-x^ - 


1 

_1 


1 

o 

1_ 

1 

8 

o 

1_ 


L 0 J 


1 

_1 


1 

_1 

1 

s 

1_ 


1 

O 

1_ 


(104) 


(105) 


(106) 
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A convention to be used to define the endpoint eonditions e will be to set 
Equation (106) equal to zero. 


e 


(x(tp) 




0 

COf 


0 


(107) 


The Mayer eost E is a form of endpoint eost. Assuming the boundary eonditions 
“must” be met, this is the measure of optimality that quantifies the solution in terms of 
the final state. Since this is a time-optimal slew, the endpoint eost is the final time. 

E[x{tf)^-tf (108) 


The Lagrange eost F is a form of running eost. This measure quantifies the eost of 
the maneuver itself A eommon form is the quadratie eost F= whieh 

represents the “energy” or effieiency of the system. Sinee the eontrol u is direetly related 
to torque, this eost ean be used to minimize torque-related funetions like power during a 
maneuver. For the simplified case of individual time-optimal slews, the Lagrange eost 
will not be used. The superseript (D.N.E.) will be used to show that the term is not only 
zero, but it does not exist for this problem. 

{^F{x(t),u(t))dt = (109) 

The cost functional 5 can be eomposed of both Mayer and Lagrange eosts. 


5[a:(-),m(')T/] = E[x{t^)^+^’ E[x{t),u{t))dt 


( 110 ) 


Rewriting this problem in the format used by Ross [25, p. 43], we get 
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Given 

Minimize 




[.] 




Subject to 


(a,,(o,,t,)^(o,o,o) 
(0/-0„m^) = (O,O) 


( 111 ) 


2. The Hamiltonian 

For every state, a costate will be defined with subscripts to match their respective 

states. 


A, 




( 112 ) 


The control Hamiltonian is given by 

ff (A,x,u) = F(x,u) + A^f(x,u) 


q(d.n.e) 

^0 ^(0 

CO 



U 


(113) 


(114) 


= Xq(0 + AjyM 


3, The Hamiltonian Minimization Condition 

The Hamiltonian Minimization Condition states that the Hamiltonian H must be 
minimized over the control u [26, p. 108]. This leads to the Euler Lagrange equation. 


du 


f=” 

(Aeft> + A^M) = A^ = 0 


(115) 


Since u does not appear, Pontryagin’s Principle must be used. From [26, 
Algorithm 3.8.13], Pontryagin’s Principle states the following. 
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dH 

du 


> 0 , 
= 0 , 
< 0 , 


—u <u <u 

max *'*'max 


U = U^ 


( 116 ) 


Applied to this problem, the switehing funetion in Equation (116) beeomes 




> 0 , 
= 0 , 
< 0 , 


-u <u<u 

max ‘'‘'max 


U = U^ 


(117) 


By rearranging Equation (117) in terms of the eostate to solve for the eontrol, the 
rule in Equation (118) results. Sinee the value of the eontrol does not matter when the 
eostate = 0, the eontrol ean be set to zero. This is not neeessary, but it is aeeeptable 
sinee the use of any eontrol during this state will not affeet the optimal eontrol problem. 
This is stated simply as an energy eonservation teehnique. Thus, the eontrol is given as: 


-Mmax \ > 0 

0 for A^=0 


Mmax < 0 


(118) 


Ultimately, Equation (118) states that the eontrol must be at its positive or 
negative maxima. The funetion sgn will be defined as the sign of its independent variable. 


sgn(.): 


1 (•)>o 

0 for (•) = 0 

-1 (.)<0 


(119) 


Equation (118) ean now be written eompaetly as 

M = -sgn(A„)M^^ 


( 120 ) 


Equation (120) matehes the formulation by [25, Algorithm 3.36]. 


4, The Costate Dynamics and the Adjoint Equation 

Erom [25, Algorithm 2.2], the adjoint equation states 
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( 121 ) 


dx 


( 122 ) 


The differential eostate equations shown in Equation (122) can be integrated to 
create equations for the costates in terms of unknown constants. 


-Cjt + C2 


(123) 


Transversality 


The endpoint Lagrangian is given by [25, Algorithm 2.9]. 


E(x^. (A:(t^ )) + v'^e[x{t ^-)) 


(124) 


E{Xf,tf) = tf + 










(125) 


The terminal transversality conditions are given by [25, Algorithm 2.8]. 




(126) 


A(a) = 


(127) 


or more simply A. = V 


(128) 
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Since the final conditions of the costates simply equal new unknown variables, 
nothing new is learned by analyzing the transversality condition. Since the final time is 
free, the Hamiltonian Value Condition is necessary. This states terminal conditions for 
the Hamiltonian [25, Algorithm 2.10]. 



(129) 


(130) 

The Hamiltonian is minimized when the dH ldu = 0 

and d^Hldu^>0. 

Equation (130) provides a fixed point for the terminal point of the Hamiltonian. Since 

dH{t) ldu = 0 , then //(r) = -1 for all time. 

6, Bang-Bang Control 


Ultimately, the primary lesson learned above is 


M = -sgn(A„)M^^ 

(131) 

K=-Cit + c^ 

(132) 

Equations (131) and (132) can be combined as 


M = -Sgn(-Cit + C2)M„^. 

(133) 


Since is linear, it ean only intersect = 0 once. This may or may not oecur 

within the range [o,t^], therefore it can be stated that at most, only one sign change 

(switch) occurs within (133). [19, pp. 193-200] additionally confirms that only one sign 
ehange occurs for initial rest-to-rest maneuvers. For arbitrary initial or final eonditions, it 
is possible for no switehes to oecur (e.g., high initial velocity with constant deceleration 
to slow down). Therefore, for the rest-to-rest maneuver in this example, the quickest 
method of rotation about this single eigenaxis would be to use bang-bang control where 
one switch occurs between [o,r^]. This means to use maximum torque available along 
that axis until half way there, then reverse torque in order to slow down to zero right at 
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the point in which the desired slew angle 6^ is achieved. The resulting trajectory is 
shown in Figure 26. 



Figure 26. Acceleration Limited Eigenaxis Slew 


Since spacecraft attitude reorientations may also be momentum-limited in 
addition to torque-limited, large angle maneuvers may look closer to that in Figure 29 
(heritage maneuver from Chapter III). Momentum is proportional to angular rate, 
therefore the trajectory is shown as rate-limited. The maneuver is initially torque-limited 
'^max(^max) ^^^^1 thc maximum momentum is reached. During this 

momentum-limited section, the spacecraft coasts at , until it can decelerate at its 
maximum torque . Figure 27 is the solution to the optimal control problem 

with an additional rate limit. 
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Torque/ Momentum Limited Maneuver 






Figure 27. Rate Limited Eigenaxis Slew 


B. OPTIMAL CONTROL EQUATIONS OF MOTION 

Since the input is already going to be determined offline per the shaped heritage 
maneuver, the time-optimal solution should instead be utilized. With theoretically 
unlimited computing power at the ground station, the optimal solution can be computed 
for every individual slew. The equations of motion for the spacecraft and reaction wheel 
system from Chapters II and IV will now be consolidated. Although the equations are 
listed in matrix form for simplicity, the number of dynamic equations reflects the number 
of independent equations that can be extracted from the equation set. 


1, Spacecraft Kinematics 

The spacecraft kinematics are repeated from Chapter II. Since the orbital motion 
is set to zero (^ = 0, see Chapter III), then . This simplifies the dynamics 

slightly allowing the removal of three states. This change is reflected immediately in 
Equation (134). 




( 134 ) 


where 


123 


^1 ^2 ^3 
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2 . 


Reaction Wheel Dynamics 


The torques of the individual reaetion wheels follow from Equation (71). 

‘'w tv*** 

Equation (135) ean be solved in terms of the reaetion wheel aeeeleration. 

*** w w 


( 135 ) 


(136) 


3. Spacecraft Dynamics 

The spaeeeraft dynamies were developed in Chapters II and IV. 


T^req=Js/c (0 + (0 CD 


N ^ W N B ^ 

req 3 w ^ 3 w 


Combining Equations (137) and (138) yields the following. 

T N -B , j N^B b rjw N^B..br^wj 

m + m xj , CO =- Z,t - co x Z.J£2 


j N -B bryw^W J N ^^B , brjwj 

Jsic ® ® X [Jsic (O + J„Q j 


(137) 

(138) 

(139) 

(140) 


Solving Equation (140) in terms of the spaeeeraft angular veloeity rate yields 


"eo' 






N ^^B . b ryw J /~\W \\ 

(O + Z,J^£2 )j 


(141) 


4, State and Control Variables 

The state veetor now has 7+n states and the eontrol veetor has n states, where n 
refers to the number of reaetion wheels in the array. 
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^1 

^2 

^3 


^4 

N ..B 

(0. 


N..B 

COy 

N^^B 

W 


n 


and u: 


(142) 


n 


n 


5, Actuator Constraints 


Torque and momentum limitations will be included in the optimal control 
problem, as these constraints represent the physical limitations of the attitude control 
system. 


T ST- ST 

max i max 


i = 1,2,...,n 


(143) 


-Q 


max i max ’ 


i = 1,2,...,n 


(144) 


6. Initial and Final Conditions 


The spacecraft will be conducting rest-to-rest slews. This means initial and final 
quaternions are known. Initial and final spacecraft rotation rates are zero. For this 
analysis, initial and final reaction wheel rates will be considered zero. In reality, 
engineers operate the reaction wheel at some nominally non-zero value. This is because 
=0 can cause jitter in the presence of “stiction” (friction from static positions) [13, 
p. 5]. The initial and final states can be represented as follows. 




N..B N..B N 

COx , (O^ 

•^0 To 




[% ’■ ■ ■ " (^ 1 ^ ’^2 ’^3 M ^ 0 , 0 , 0 , 0 ,- • ■ , 0 ) 


(145) 

(146) 


65 



7. 


Cost Function 


The objective is to minimize the overall maneuver time. Therefore, the cost 
function is 

5[jc(-),M0T/] = ty (147) 


8, Optimal Control Problem 

By gathering Equations (134) to (147), the optimal control problem can now be 
defined; 


Given 


x = 


N 

q, ^3 q^ (O^ 


CO. 


(0^ Qi 


Q 




Minimize 





Subject to ^123 = ^[q, ^(o’" - ^(o’" x 

1 / JV B 

^4=-2l ® 







W N B 

W _ N H 

s/c 


+ ''z;7„i2"')) 





a'= j:'< 





{% 




, CO^ , co^ . 

^0 To ■ 

^^^B r^W 

r\W \ / 0 0 0 

,^4,0,0,0,0,---,0) 

{% 

’^2/ 

’^3/ 



N^^B r^w 

'.‘ 4 :, )=(<?/ 

f f f 

.^2 ^3 ^4 

,o,o,o,o,---,o) 





max 

<t”' <t'^ 

1 max 









max 

<a'^ <a'^ 

~ 1 ~ max 






C. SCALING THE PROBLEM 


(148) 


The purpose of scaling is to generalize all states to within similar orders of 
magnitude. Without scaling, one state may take on large values (lO^) while others are 

comparatively small (lO *’). The relative difference would be is (lO which causes 
numerical problems due to operating near machine epsilon. By simply multiplying each 
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state by the appropriate “sealing factor” their values can be made the same order of 
magnitude, and numerical issues successfully sidestepped. 

The common formulation for scaling units is done through “canonical units.” 
Canonical refers to a specific definition where, for example, a velocity unit (VU) is equal 
to the quotient of a distance unit (DU) and a time unit (TU). 

VU = — (149) 

TU 

Within [21], “designer units” are chosen instead to keep the states and costates in 
relative orders of magnitude. Instead of looking at the costates, nor at the coupling 
between position and velocity, this thesis outlines a different kind of designer unit. They 
will be denoted by a superscript * (or state followed by “z” within any listed code). 


1, Designer Units 

As scaling factors, the following designer units are used to bring all variables 
close to the first order of magnitude. 

Wheel torque: T* = (maximum wheel torque for chosen RW) 

Time: T* = Is (use real time) 

Quaternion: No further scaling. Assumes -1 < < 1 is sufficiently scaled 

S/C Spin Rate: co* = 3° I s = 0.0524 rad / 5 (slightly faster than conventional rate) 

Wheel Rate: Q* = (maximum wheel spin rate for chosen RW) 

For example, the relationship between the torque T and the scaled torque t will 
be the proportionality constant T*. 

t = Tt* (150) 

Therefore, the scaled torque would always operate in the range [-l,l]. 



( 151 ) 
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Scaling of constants within a canonical unit representation ean be useful beeause 
many terms will eaneel out leaving simpler equations to manipulate and eode. If all 
sealing terms are direetly eoded into the software, sealing of eonstants does nothing. This 
is beeause the eonstant must be scaled down first, then sealed baek up to its original value 
within the dynamie equation. In the case of designer units where many mismatched units 
are being manipulated, it is unlikely that eonvenient units will cancel leaving an elegant 
equation. Therefore, sealing is not required and not eompleted for constants. An example 
of designer unit scaling is as follows for the spaeeeraft dynamies. The full set of scaled 
equations is shown in Appendix A. 






hrjW,W 
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^ (O*^ 
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(152) 
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2. Scaled Optimal Control Problem 

The sealed optimal eontrol problem is 
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Given 


X = 


u = 


— — — — A'—A/— a^;^d 

Ql qi ^3 ^4 ^1 ■•• 

-\T 


T," 




Minimize 5 [x (•),«(•),?/]=?/ 

Subject to 9 j23 = ^(^4 X 9 j 23 )(r*(u*) 


^4 = -^(Wf9,23(rV) 
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® =-/./c 


f T*r* ^ 
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-1 < < 1 


(153) 


-i<q!^ <1 


D, DIDO IMPLEMENTATION 

DIDO is a MATLAB-based software package for solving optimal control 
problems [21], A huge advantage of DIDO is that the optimal control problem from 
Equation (153) is the only necessary derivation. Since Pontryagin’s Principle, discussed 
in Section A, is automatically integrated by the DIDO software, further analysis is not 
required. But, Pontryagin’s Principle can be checked to verify the numerical solution. 
The code simply requires the optimal control problem to be inputted into a series of four 
functions run from a single m-file. 

1. Cost Function 

The cost function comes from Appendix A, Equation (275), implemented in 
Eigure 28. 

5[^(-),«(')T/] = ^ (154) 
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% Cost function for the Min Time Slew using RWs 

function [endpointCost, runningCost] = MinTimeSlewRWCost(primal) 

tf = primal.nodes(end); 

endpointCost = tf; 
runningCost = 0; 


Figure 28. Cost Function Implemented as a DIDO Function 


2, Endpoint Function 

The endpoint function contains the boundary conditions from Equations (280) and 
(281), implemented in Figure 29. 


1 77 77 77 77 N ttB N ttB N ttB \ _ 

(— — — — N —B N —B N —B iF.W iF.W \ 


y 0 0 0 0 ^ 

^^^^0000 •••0 

* 9 * 9 * 9 * 9 '^ 9 '^ 9 '^ 9 '^ 9 9 '^ 

yq q q q j 

* 9 * 9 * 9 * 9 '^ 9 '^ 9 '^ 9'^9 9 '^ 

\q q q q ) 
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% Events function for the Minimum Time Slew using RWs 
function endpointFunction = MinTimeSlewRWEndpoint(primal) 

global CONSTANTS 

Zn = CONSTANTS .Zn ;% number of RWs 

qO = primal.states(1:3,1) ; qf = primal.states(1:3,end) 

q40 = primal.states(4,1) ; q4f = primal.states(4,end) 

wO = primal.states(5:7,1) ; wf = primal.states(5:7,end) 

WO = primal.states(8:7+Zn,1); Wf = primal,states(8;7+Zn,end) 

endpointFunction = [qO; q40; wO; WO; qf; q4f; wf; Wf]; 


Figure 29. Endpoint Eunction Implemented as a DIDO Eunction 


3, Path Function 

The path is not necessary for this problem, so no path function script is needed. 

4, Dynamics Function 

The dynamic equations are from Equations (276) through (279). 
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Within the dynamics function shown in Figure 30, state scaling occurs at the point 
of state assignment and is then accounted for in the equations of motion. 
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% Dynamics function for the Minimum Time Slew using RWs 
function xdot = MinTimeSlewRWDynamics(primal) 


global CONSTANTS 
Js = CONSTANTS. Js 

Jw = CONSTANTS .Jw 

Z = CONSTANTS .Z 

Zn = CONSTANTS. Zn 


; % spacecraft inertia matrix 
; % wheel inertia matrix 
;% RW orientation 
; % number of RWs 


Tauz 

= CONSTANTS. Tauz 

Wz 

= CONSTANTS. Wz 

wz 

= CONSTANTS. wz 

Tz 

= CONSTANTS. Tz 

qz 

= CONSTANTS. qz 


;% torque unit 
;% wheel rate unit 
;% wheel rate unit 
; % time unit 
;% quaternion unit 


% States and controls 


q 

= primal.states(1:3,:)*qz 

;% 

quaternions 1-3 

q4 

= primal.states(4,:)*qz 

;% 

quaternion 4 

w 

= primal.states(5:7,:)*wz 

;% 

spacecraft rate 

W 

= primal.states(8:7+Zn,:)*Wz 

;% 

wheel rates 

tau 

= primal.controls*Tauz 

;% 

wheel torques 


N = length(q4); %determine size and preallocate matrix sizes 
qdot = zeros(4,N); q4dot = zeros(1,N) ; 
wdot = zeros(3,N); Wdot = zeros(Zn,N); 


% Equations of Motion: 
for i = 1:N 


qdot(:,i) = l/2*( q4(i)*w(:,i) - cross(w(:,i),q( : ,i) ) )*Tz/qz ; 

q4dot(:,i)= -l/2*w(:,i)'*q(:,i)*Tz/qz ; 

wdot(:,i) = Js\eye(3)*(-Z*tau(:,i) ... 

-cross(w(:,i),Js*w(:,i)+Z*Jw*W(;,i)) )*Tz/wz; 

Wdot(:,i) = Jw\eye(Zn)*tau(:,i) *Tz/Wz; 

end 

xdot = [qdot ; q4dot; wdot ; Wdot]; 


Figure 30. Dynamics Function Implemented as a DIDO Function 


5. Main Function 

With only a few dozen lines of code, the main m-file runs the four functions and 
DIDO. This is where the CONSTANTS are defined. Boundaries are set up for the states, 
controls, and endpoint conditions (see Figure 31). 
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tfMax = 40 


;%[sec] upper bound for scaled problem 


bounds.lower.time = [0; 0]; 

bounds.upper.time = [0; tfMax); 

% Scaled initial conditions 
gO = q_0/qz ; 

wO = [0;0;0)/wz ; 

WO = zeros(Zn,1)/Wz ; 

hO = Jw*W0 ; 


% final conditions 
qf = q_f/qz 
wf = [0;0;0]/qz 
Wf = zeros(Zn,1)/Wz 

% boundary conditions [ ql-q4 ; wx;wy;wz 

bounds.upper.states = [ones(4,1)/qz; ones(3,1) 
bounds.lower.states = -bounds.upper.states; 

bounds.upper.controls = ones(Zn,1)*tau_max/Tauz; 
bounds.lower.controls = -bounds.upper.controls; 

bounds.upper.events = (qO; wO; WO; qf; wf; Wf); %(20xl) 
bounds.lower.events = bounds.upper.events; 


Wl-Wn ] 

ones(Zn,1)*W_max/Wz j; 


% Tell DIDO the bounds on the problem 
MinTimeSlewRW.bounds = bounds; 


Figure 31. DIDO Bounds Implemented within a DIDO seript 


DIDO is initially run with no guess and as little as 20 nodes to determine if a 
feasible solution will arise (see Figure 32). 


%% RUN 1 - NO GUESS MADE 

% Number of nodes for the spectral algorithm 
algorithm.nodes = 20; 

[cost, primal, dual] = dido(MinTimeSlewRW, algorithm); 

Figure 32. DIDO No Guess (20 nodes) Implemented within a DIDO script 

The results are then seeded back into DIDO as initial guesses along the entire 
dynamic range using the following code. This time a much more refined set of 120 nodes 
is used (see Figure 35) 
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%% RUN 2 - SEED RESULTS OF RUN 1 AS GUESS 

algorithm.nodes = 120; 

algorithm.guess = primal; 

[cost, primal, dual] = dido(MinTimeSlewRW, algorithm); 

Figure 33. DIDO Seeded Guess (120 nodes) Implemented within a DIDO script 

This completes all steps required to run DIDO under both the “no-guess” and 
“seeded guess” scenarios. The results from run 1 can be looked at as an initial estimate to 
put DIDO on track. 

E. OPTIMAL MANEUVER IMPLEMENTATION 

The optimal control solution is fundamentally open loop and the control solution 
is applied directly to the individual wheels (see Figure 34). Additionally, the optimal 
quaternion path is presented as a shaped angular trajectory at the forefront of the loop. 
Since body rates are known, the EQF can follow the rate trajectory instead of being 
regulated to zero angular velocity. When plant uncertainties or external disturbances 
exist, the instantaneous quaternion error generates work for the EQF Controller, applying 
any additional wheel torque as required. With these changes. Figures 2 and 17 become 
Figures 34 and 35, respectively. 


OMal Mobon 



Figure 34. Architecture for Optimal Maneuver Implementation: Spacecraft Model 
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Figure 35. Architecture for Optimal Maneuver Implementation; Reaction Wheel 

Subsystem 


Figure 36 shows the integration of the rate error into the eigenaxis quaternion 
feedback controller, a slight modification from Equation (28). 

^req ~ Js/c9ei22~ ^s/c ® ^ Jsic ® (157) 



Figure 36. Architecture for Optimal Maneuver Implementation: Eigenaxis Quaternion 
Feedback Controller 
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By modifying the control architectures as per Figures 34 to 36, the full reaction 
wheel array torque and momentum is available for use. If the state quaternions q feed 
back disturbances and/or uncertainty error, quaternion error is corrected by creating the 
quaternion error from the optimal quaternion trajectory. The EQF Controller then creates 
a torque requirement. Body rate error is propagated through both the EQF controller and 
wheel dynamics in the form of additional required torque. The total required torque is 

then passed through the pseudo-inverse to sum with the optimal torque . 


X 


ir 

W 



(158) 


The system will naturally require some correction torque xf ^^, because the 

spacecraft inertia and external forces (if present) cannot be modeled perfectly. This leads 
to the requirement for the pseudo-inverse to perform the allocation of the feedback 
torques. Since the pseudo-inverse uses the minimum least-squares solution for converting 
body torque into n wheel torques, the feedback provided follows a minimum energy 
control allocation. 


To accommodate the need for feedback, it is necessary to reduce the control 
authority of the optimal control solution. This is because if the optimal torque solution 

is already utilizing the maximum system torque x '^^, the system will not have the 

authority needed to make the feedback corrections. The available control authority for 
optimal control should be reduced inversely proportional to the knowledge of the system. 
If the system knowledge is high, then the control authority should be reduced only 
slightly. For instance, if the system knowledge is high, an optimal control authority of 
95% may be appropriate. If the system knowledge is low, an optimal control authority of 
75% may be appropriate. Within the optimal control problem formulation presented here, 
the adjustments would resemble that in Equation (159). 


-Q.95xZ^<x'f <0.95 X 


W 

max 


(159) 
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F. EXAMPLE TIME-OPTIMAL MANEUVER 

In this section, the feasibility and optimality of a representative time-optimal 
example will be documented. As a representative generic spacecraft, the following 
figures will be used. 

-w 

0 -V2 0 

V2 0 -^/2 (160) 

1 1 1 


2.54 0 0 T 


0 2.54 0 kg-rn^ 

(161) 

0 0 2.54 


^u,,max = ±8.57mA-m 

(162) 

= ± 3000 rpm = ±314.2raJ/5' 

(163) 



A 0 0 0 

•/W, 33 

0 A 0 0 

J Wj 33 

0 0 A 0 

0 0 0 j„ 

VV4 33 

The designer units selected are: 

Wheel torque: T* = ±8.57mA-m 

Time: r* = lsec 



(164) 


Quaternion: ^* = 1 (no further scaling) 

S/C Spin Rate: co* = 3deg/ s = 0.0524rad / s 

Wheel Rate: Q* = 314.2raJ / s 
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A 120° rest-to-rest slew will be performed about the eigenaxis 


e = 


^[ 2/3 0 VTTs 


from q° = 


-Vr72 0 -1/2 1/2 


-\T 


to the origin 


9^=[o 0 0 l] - The boundary conditions are shown in Equation (165). 

— N—B N—B N—B \ 

= (-Vr72, 0,-l/2, l/2,0,0,0,0,0,0,0) 


(165) 


:( 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 ) 


1, Optimal Solution 

The time history of the control allocation determined by DIDO is shown in Figure 
37. For each of the 120 time nodes, DIDO determined the optimal control solution for 
each wheel. Finear interpolation was used to determine the feedforward control allocation 
that would be used for the simulation. 


RW Torque (mN-m-s) for [-0.70711,0,-0.5,0.5] to [0,0,0,1 ] slew (120 nodes) 
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Figure 38 shows the time-optimal quaternion trajeetory. The final node exists at 
= [ 0 0 0 1 ] as expeeted in the boundary eonditions. 


Quaternions for [-0.70711,0,-0.5,0.5] to [0,0,0,1] slew (120 nodes) 



Figure 38. Optimal Quaternion Trajectory (120 Nodes) 

Additionally, the final nodes for the angular rates all end at zero per the boundary 
conditions. 

2, Feasibility 

A feasible control solution is one that propagates through the system dynamics as 
expected and satisfies the boundary conditions. The feasibility of both solutions is 
checked by propagating DIDO’s control solution through the system dynamics. First, the 
state values at each of the individual time nodes must align with the propagation of the 
dynamics at those specific times. Second, the dynamics must propagate to the appropriate 
final conditions. If both of these have been satisfied, then it is determined that enough 
nodes have been used and the solution is indeed feasible. 
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The interpolated control signal (Figure 37) is propagated through the spacecraft 
dynamics, modeled within Simulink per Figures 34 to 36. For DIDO propagation, EQF is 
turned off. The DIDO torque signal is fed into the reaction wheel dynamics for open-loop 
propagation. Since the knowledge of the spacecraft inertia is considered perfect, the 
feedback portion will only affect the conventional slew, later used for comparison. 

States: Once propagation is complete for the full 120 nodes, the control and states 
for both runs and their respective control propagations can be plotted together. If the 
curves for each propagated state pass through its respective DIDO nodes, the DIDO- 
produced control solution can be deemed feasible. 

Figures 39-41 show is a graphical comparison of the initial DIDO solution vs. the 
dynamically propagated solution from Simulink for all 11 states. By inspection, it appears 
that every single numerically calculated point is directly on top of its propagated 
trajectory. This shows that the DIDO solution is not only feasible, but also very accurate. 


Quaternions for [-0.70711,0,-0.5,0.5] to [0,0,0,1] slew (120 nodes) 
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(deg/s) for [-0.70711,0,-0.5,0.5] to [0,0,0,1] slew (120 nodes) 



RW Rate (rpm) for [-0.70711,0,-0.5,0.5] to [0,0,0,1 ] slew (120 nodes) 
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Scaled States and Costates: Another important check to determine the efficiency 
of code is the relative magnitude of the scaled states and costates. In Figures 42 and 43, 
the scaled states and costates {x,X) are represented in the selected designer units prior to 
their rescaling back into SI units. All states and controls are on the same order of 
magnitude, which demonstrates good scaling. Only a 1 order of magnitude difference 
shows up in the costates, which are very sensitive to scaling. The relative comparison 
shows that reasonable scaling was achieved during problem formulation. If the costates 
were to differ from each other by 1-2 orders of magnitude, scaling of the time constant 
T* or turning the knobs on the state parameters could adjust this. Additionally, a key 
point to note is that the scaled states only differ from their costates by one to two orders 
of magnitude. A well-scaled problem will have them roughly the same order of 
magnitude [21, p. 33], therefore tweaking of the designer units could remedy this. 

1 

D" 0 

-1 

1 

3 0 

-1 

1 

Ci 0 

-1 

0 5 10 15 20 25 30 35 40 

time (sec) 

Figure 42. DIDO Scaled States 
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Scaled Costates for [-0.70711,0,-0.5,0.5] to [0,0,0,1] slew 







Figure 43. Scaled Costates 


3, Optimality 

The first place to check for optimality is the Hamiltonian. To ensure the optimal 
solution is found, the Hamiltonian must be held constant with respect to time. 
Additionally, the final time is free, so the Hamiltonian Value Condition is necessary, 
refer back to Equations (124) through (130) from Section A. Similarly, it can be shown 
that 

= = ( 166 ) 

Equation (166) reveals a key point; for all time-optimal cases, the Hamiltonian is 
held constant at -1. Eigure 44 is a zoomed-in plot that shows that the time-history of the 
Hamiltonian on approximately -1 ± .007, which is reasonable. 
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Hamiltonian 



Figure 44. Hamiltonian 


Although DIDO does not require inputs from the optimality analysis in Seetion A 
to run, this information ean help further analyze the optimality of the solution. 
Pontryagin’s switehing funetions are easily derived. For the four eontrol torques, the 
switehing funetions are as follows. 

5, = 4, / A. - VTTI4, / A - ^/273 /1 
s, = 4, / 4 - VTTI4, / A - ^/273 / A 

I — /- 

4 = 43 / 4 - VI73 a„3 / J, - 7573 4, / j, 
s, = 4. / 4 - 7173 4, / A+7773 / /, 

where 

Per Pontryagin’s Prineiple (refer to Seetion A), these switehing funetions 
represent the slope of the Hamiltonian with respeet to the eontrol used. This defines the 
eontrol spaee by further determining what eontrol effort should be used at any instant in 
time. If the slope is negative, i.e., S.<0 , the maximum torque should be used (minimum 

torque for positive slope). When the slope = 0, the value of the eontrol does not really 
affeet the Hamiltonian. The following switehing logie represents the optimal use of the 
eontrols. 


If sA 


<0 ~ '^max 

= 0 , then use T„,„<T,<T„ 
>0 T:= 


(169) 
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Figure 45 shows the dynamics of the reaction wheel assembly alongside the 
solution for the switching function at each of the 120 nodes. The torque responds very 
closely to the switching function logic in Equation (169), thereby further validating the 
optimality of the solution. 


RW Torque (mN-m-s) and Switching Function (120 nodes) 



G. EIGENAXIS SLEW COMPARISON 

For comparison purposes, the Eigenaxis Quaternion Feedback (EQF) Controller is 
used to perform a conventional eigenaxis slew maneuver via a shaped bang-coast-bang 
trajectory. EQF is governed by Equation (170), repeated from Chapter III, Equation (28). 

T., = -cJ(D+^(0^x (170) 

Eor this example, a 0.1 sec settling time was chosen. Recall, the step size is a 
function of the distance between discrete points used within the shaped quaternion 
trajectory. Therefore, a small settling time can now be used. The damping ratio was 
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chosen as 0.9 so that it does not allow for much overshoot. Using these criteria, the gains 
k and c can be found. 

^..«fe=0-lsec 

C = 0.9 

®„ = 4/(t,.„,0 = 44.4 (171) 

k = a)l = \.91E^ 
c = 2C«„ = 80 

A sanity check can be done by comparing the time-optimal slew against the 
classic EQF mentioned above. Figure 46 demonstrates the different paths taken by each 
slew. Notice that q 2 changed for the optimal slew only. Acceleration also occurred faster 
for the optimal slew. A detailed analysis of how acceleration and momentum are utilized 
differently occurs within Chapter VII. 



time (sec) 

Figure 46. Optimal and Figenaxis Comparison; Quaternions 

Figure 47 represents this comparison with respect to the Euler angles. The DIDO 
run took only 40.5 seconds, 21.8% faster than the 51.8 second EQF slew. It does not 
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matter which axis is used, the time-optimal slew will always be just as fast or faster than 
the eigenaxis slew. 


EQF Slew Comparison (Optimal Slew is 21.8% faster) 



Figure 47. Optimal and Eigenaxis Comparison; Euler Angles and Rates 


The optimal maneuver achieved the final boundary conditions much faster for two 
reasons. It was able to take advantage of off-eigenaxis motion, achieving much greater 
magnitudes in both torque and momentum. Eigure 48 illustrates the additional torque 
(steeper slope) and momentum (maximum rate). The body rate of the eigenaxis slew 
confined by the pseudo-inverse was restricted to 3°/s. On the other hand, the optimal 
maneuver could reach momentum and torque outside the pseudo-inverse space, obtaining 
body rates of up to 4.24°/s, a 41.5% increase in maximum rate. 
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Figure 48. Optimal and Eigenaxis Comparison; Body Rates 


Table 2 presents a maneuver-time eomparison of a few more 120° slews. Each is 
moved from the initial quaternion back to the origin [0,0,0,!]. As the initial quaternion 
changed, eigenaxis slew rate increased. This is because the eigenaxis was limited only by 
the pseudo-inverse-limited octahedron, not the spherical surface (see Eigure 24). As the 
eigenaxis approached a maximum vertex in run 2, the optimal slew only improved the 
time by 4.2%. Interestingly, the time-optimal maximum body rate was consistently 
between 4.0-4.3°/s, but this accounts for much off-eigenaxis movement. 


qo 

Eigenaxis 

Eigenaxis Slew 

time 

(sec) (7s) 

Optimal Slew 

time 

(sec) (7s) 

time 

improvement 

[-0.707,0,-0.5,-0.5] 

[0.8165,0.0000,0.5774] 

51.8 

3.0 

40.5 

4.2 

21.8% 

[-0.5,-0.5,-0.5,0.5] 

[0.5774,0.5774,0.5774] 

44.0 

3.7 

42.1 

4.3 

4.2% 

[-0.612,-0.433,-0.433,0.5] 

[ 0.7071, 0.5000, 0.5000] 

46.7 

3.5 

41.5 

4.0 

11.0% 

[0.612,0,0.612,0.5] 

[-0.7071, 0.0000,-0.7071] 

51.2 

3.1 

39.8 

4.3 

22.3% 


Table 2. Optimal vs. Eigenaxis Comparison (120° slews) 


H, EFFECTIVE EIGENAXIS OF TIME-OPTIMAL MANEUVER - 
ANALYTICAL APPROACH 

It is important to introduce the concept of effective eigenaxis. Given any two 
quaternions, the optimal slew deviates off-eigenaxis, varying its rotation axis and rate 
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throughout the maneuver. The effeetive eigenaxis is defined as the axis of rotation for a 
minimum angle slew (original eigenaxis from qo to qf), regardless of path aetually 
traveled during the maneuver. For example: for Table 2, run 1, the effeetive eigenaxis of 
the optimal slew was [0.8165,0,0.5774]. The aetual distanee traveled by the time-optimal 
slew was 123° versus the minimum angle maneuver (120°) about the eigenaxis. 

Using this concept, it is possible to determine the relationship between the 
effective eigenaxis and the effective eigenaxis torque and slew rate about that effective 
eigenaxis, the motivation for the next two sections. This is necessary to form a general 
comparison between the two very different maneuvers. The crucial advantage of the 
optimal maneuver is that it can access all of the dodecahedron momentum space, 
including that outside of the pseudo-inverse octahedron (Figure 49). This is because the 
optimal torque trajectories for each wheel are obtained as part of the solution to the 
optimal control problem. Flence the pseudo-inverse need not be employed (recall control 
law from Figure 34). This momentum envelope minimum begins at 1.633 at the 

center of the twelve flat faces, expanding to 2.309 at the maximum vertices (1.633 
to 2.309 in torque space). 



Time Optimal 
can access 
these previously- 
unavailable 
regions 




Figure 49. Time-Optimal Accessible Regions 
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It is difficult to compare a generic optimal maneuver against an eigenaxis slew, 
primarily beeause the precise shape of the optimal solution is not intuitive. The optimal 
maneuver is tailored to maximize performance of an individual maneuver. Thus, a 120° 
slew about one effeetive eigenaxis might be far different than a 120° slew about a 
different effeetive eigenaxis (see Table 2 for slew time comparison). Looking at Figure 

50, the [ 0 1 0 ] (or b^) axis yields a limit at 1.633 ('Z^max)’ while the eigenaxis 

about 1^ 0 0 1 ] (or hj) yields 2.309 ('Z^max)- This means that the rotation about 

the ^3 axis ean potentially accelerate and traverse 41% faster than a rotation about . 

Here, the effeetive eigenaxis momentum will be derived from the average 
eigenaxis momentum over 47i steradians. The first step will be to caleulate the average 
amplitude of momentum or torque available about every possible eigenaxis. The average 
amplitude of a series of any set can be determined by the sum of the set divided by the 
number of entries in the set. Therefore, the average of magnitude of n (100, 1000, IM) 
eigenaxis momentum vectors is the quotient of the sum and n. 



If was placed back into 3D space along the same n eigenvectors, the shape 

would now be a sphere of radius . Consider eaeh veetor to have some volume 

assoeiated with it. Since the total length of the sum is eonserved, the sphere would have 
the same volume as the original dodecahedron. 


This approach will be followed in reverse in order to attain , the radius of the 

sphere. Using MATLAB’s convhull eommand, the volume of a shape is caleulated as a 
subset of its three-dimensional points. Using this command, the volume of the 
dodeeahedron is caleulated as 7i4.632 . A simple spherical relationship yields the 

radius of the sphere. 



( 173 ) 
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Now, the equivalent spherieal radius of the reaetion wheel momentum (or torque) 
envelope is 1.805. The shape is shown in Figure 50. 


h , , =1.805/1 

opt-sphere rr 

'^opt-sphere ~ ^ -805 


(174) 


Since the sphere is not inscribed within the dodecahedron, it would appear that 
some torque/momentum described by the sphere is not available. This may not be the 
case though as the sphere simply represents the effective torque/momentum about the 
effective eigenaxis. Since the trajectory is allowed to deviate from the effective eigenaxis, 
the torque/momentum space at the vertices is exploited to raise the effective 
acceleration/rate limits outside the boundaries of the torque/momentum envelope. 
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Consider the optimal trajectory using only eigenaxis slews, but has access to the 
entire momentum space. The average slew rate (or acceleration) derives directly from the 
newly-created optimal momentum sphere with values per (174). Since both torque and 
momentum limitations are assumed, the average time for any given slew angle can be 
predicted using the simple trapezoidal trajectory: bang-bang (from Figure 26) or bang- 
coast-bang (from Figure 27), depending on whether the momentum limit is reached. The 
same spacecraft/wheel configuration from Section G will be used for comparison. Over 
the slew regime extending to 180°, the curve shown in Figure 51 is developed. 



Figure 51. Reaction Wheel Effective Eigenaxis for Time-Optimal Maneuver - Analytical 

0s=25.4 kg-m^) 


The curve at the beginning indicates a torque limited region. Eor slew angles of 
30° or less, the reaction wheel array does not have the torque to reach its full momentum. 
Therefore, the maneuver is bang-bang. It will be in maximum torque mode (edge of the 
torque sphere, 0.1805Am), then reverse maximum torque to slow down without ever 
reaching the edge of the momentum envelope. Eor slews greater than 30°, the maximum 
momentum is reached and the system theoretically coasts at a constant rate (proportional 
to 1.805 Nms ) for a portion of the maneuver. 
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In this section, the concept of effective eigenaxis was used as a baseline to 
compare the full spectrum of possible optimal maneuvers against the predictability of 
shaped eigenaxis slews. The relationship between effective eigenaxis torque/momentum 
for optimal maneuvers was developed. For the optimal maneuver, an effective 
momentum of 1.805 and torque of 1.805 can be used to represent the potential of 
the time-optimal maneuvers over 47i steradians. Additional comparisons will be carried 
out in Chapter VII. 


I, EFFECTIVE EIGENAXIS OF TIME-OPTIMAL MANEUVER - 
NUMERICAL VERIFICATION OF ANALYTICAL APPROACH 


The analytical approach from Section H must be validated. To gather sample 
data, 100 optimal simulations were completed using the spacecraft/wheel configuration 
from Section G. The only dependent variable used was the beginning attitude. A random 
quaternion was generated for each trial, and the system slewed to the origin. 


= 2(ran(i(4,l)-0.5) 

= [ 0 0 0 1 7 


(175) 


This allows a random sampling of slew angles and eigenaxes. The sample data is 
plotted alongside the analytical curve, shown in Figure 52. 
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Figure 52. Reaction Wheel Effective Eigenaxis for Time-Optimal Maneuver - Numerical 

Approach with E807 Xmax 0s^25.4 kg-m^) 


The majority of the sample data lies above the curve, demonstrating that the 
approach taken in Section H might be too aggressive. The slope of the momentum- 
limited region seems correct, but the intercept with the torque-limited region is too early. 
Therefore, the effective torque bound must be reduced. Several values for the effective 
torque are shown in Eigure 53. 
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Figure 53. Reaction Wheel Effective Eigenaxis for Time-Optimal Maneuver - Numerical 

Approach with [E333,E569,L807] Xmax Gs=25.4 kg-m^) 

The sample data points in Eigure 53 are seen to lie below the curve, 

demonstrating that the spherical system devised during the analytical approach made an 
achievable configuration. Since the actual optimal runs were faster than the curve fit, it 
shows that the spherical approach using the reduced value of tends to give a slightly 

conservative estimate of the true performance of time-optimal solutions. This 
demonstrates that the maneuver times obtained by assuming a trapezoidal input via an 
effective eigenaxis slew is achievable via optimal control all the time. This allows the 
spherical shape with momentum/torque radii equal to E805/E333 of the individual wheel 
momentum/torque maximum to be used in later in this thesis to facilitate comparison. 
This 35% increase in momentum over 47i steradians is significant. 
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opffl, 


1.805/z„ 

1.333/z„ 


= 1.354 


(176) 


In Chapter Vll, it will additionally be shown how the optimal maneuver 
maximizes the use of the momentum/torque envelope to aehieve better-than-eigenaxis 
performanee. 


J. EFFECTIVE EIGENAXIS OF PSEUDO-INVERSE LIMITED TIME- 
OPTIMAL MANEUVER 

If adding optimal torque signals directly to the individual wheels is out of the 
question, then the system can still be optimized for trajectory shaping by using a slight 
variation from the time-optimal logic used in Figures 34-36. The system must still 
employ the pseudo-inverse within the reaction wheel dynamics, but performance can be 
improved within this space. Figure 54 illustrates the pseudo-inverse limited optimal 
shaped trajectory being used in the spacecraft model. 
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Figure 54. Pseudo-Inverse Limited Time-Optimal Spacecraft Model 


The accessible gains are highlighted in Figure 55. Granted, the blue regions are 
accessible without the use of optimal control, but not about every eigenaxis. This section 
demonstrates how optimal control can be exploited for equal gain over 47i steradians. 
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Figure 55. Pseudo-Inverse Limited Time-Optimal Accessible Regions 


In the same fashion as analytically extracted for the time-optimal in Sections H 
and I, the analytical and numerical approach can be completed for the pseudo-inverse 
limited time-optimal maneuvers. It can be shown that the pseudo-inverse limited time- 
optimal space can be achieved by using maximum magnitudes of 1.577/1.333 of the 
individual wheel momentum/torque maximums. 


h =1 511 h 

pseopt-sphere * max 

^pseopt-sphere 1 *333 ^ 


(177) 
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Figure 56. Pseudo-Inverse Limited Time-Optimal Momentum Sphere 


By shaping the trajectory using pseudo-optimal inputs rather than the standard 
trapezoidal input, 18% increases in speed and acceleration can be achieved. 


pseopt^^ 


B 


\.511h„ 

I333h 


1.183 


(178) 


K. SUMMARY 

This chapter was used to develop and analyze the time-optimal slews. First, the 
time-optimal slew for the eigenaxis-restricted case was demonstrated as a bang-bang, or 
as a bang-coast-bang maneuver when momentum-limited. Next, the eigenaxis was 
abandoned and the time-optimal reorientation was developed and implemented in DIDO. 
Slew time improved over the trapezoidal input by up to 22% for 120° slews. This optimal 
control approach allowed the limitations of the pseudo-inverse to be avoided, increasing 
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the available momentum and torque envelope. Next, the effeetive eigenaxis concept was 
introduced, which allowed a comparison between the eigenaxis reorientation and the 
time-optimal slew. The effective eigenaxis was further developed to show how, on 
average, the available momentum and torque for the time-optimal case are 1.805 and 
1.333 about the effective eigenaxis. For the pseudo-inverse-limited time-optimal 
case, they are 1.577 and 1.333 , respectively. 
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VI. CONTROL MOMENT GYROSCOPES 


Like reaction wheels, CMGs also have their physical limitations. They have 
maximum angular rates and accelerations that must be modeled. This ensures that only 
the maximum individual CMG torques can be applied, regardless of demanded torque. In 
order to appropriately apply limitations, a model of the CMGs must be created. This 
chapter is a focus on single gimbal control moment gyroscopes (SGCMGs). Dual gimbal 
and variable speed CMGs will not be discussed. The following introduction is simplistic 
regarding SGCMGs, but representative of the key interactions. The rest of this chapter 
will be used to derive these terms and show the interaction with the spacecraft body. 

At the center of a CMG is a large momentum wheel, driven by its own wheel 
motor at a nominally constant speed. The basic concept of a CMG is for a wheel to spin 
freely within a secondary gimbal frame. The gimbal frame is rotated within the body 
frame by a gimbal motor. The change in the direction of the wheel angular momentum as 
a result of gimbal motion causes a large torque orthogonal to the gimbal and angular 

momentum axes, i.e. f = <5 x h . Since a relatively small gimbal torque input can become 
a much larger output torque, this concept is known as “torque multiplication” [27, p. 1] or 
“torque amplification” [28, p. 160]. 

A basic model of a CMG is shown in Figure 57. At the center of the CMG is a 
momentum wheel shown in the wheel frame W. The Wj and H ’2 are fixed to the 
wheel and rotate in the gimbal frame as a function of the constant wheel rate Q and time. 
It is important to note that the spin axis of the wheel frame is aligned with the gimbal 
frame, such that 

=Q.w, = ag, . (179) 

The wheel motor is attached to the gimbal frame and drives rotation of the wheel. 
The wheel is spun at a constant rate; therefore, the simple rotation of the wheel alone 
does not cause a torque transfer between the wheel and gimbal frames. 

< = 7„Qg3=0 (180) 
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Figure 57. CMG Schematic 


Figure 57 shows the small gimbal momentum and torques as small red and green 
arrows, respectively. The large wheel momentum and large resultant orthogonal torque 
are shown in with larger arrows. All of these momentums and torques can be labeled 
relative to the gimbal frame. For simplicity, gyroscopic coupling terms are not shown on 
the diagram. 


K = Ksi 

'^w ~ '^wSl 


(181) 

(182) 


Consider the spacecraft body frame aligned with the gimbal frame shown in 
Figure 57. Since the gimbal frame is allowed to rotate within the body frame, the 
following set of time derivatives with respect to the body frame are important to note. 
They will be used in the chapter’s derivations. 
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(183) 


dt 


(^l) = <5g2 


dt 




dt 


(^ 3 ) = 0 


A. ORIENTATION 

The description of the orientation of a CMG will be carried out much like that of 
the orientation of reaction wheels in Chapter IV. To avoid redundancy, some portions of 
the explanation will be condensed. 

The fundamental formula for the angular momentum gimbal to body change of 
basis is given as 


The dyadic form of this projection onto the body frame is 


K = \bgl]hr‘ = 


k*k k*k 


h! 


(185) 


The orientation will be placed in the B frame, converting to the G,. frame via 
, and then using the transpose to convert back. In the example shown in Figure 
58a, the MED is placed directly on the B xyz origin, rotating about the axis and 
gimbaling around the k per Figure 58b. This produces an angular momentum 
directly in the k direction for the wheel. 
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a) Initial aligned 
with body axis 



b) After 3 (a) ^ 3 a ^3 

rotation 


k 


^2 ^2 


c) After 2 (P) 
rotation 



d) After 3 (6) 
rotation 



Figure 58. CMG Transformation 

The first two orientation transformations are exactly like those for the reaction 
wheel. A fundamental 3 rotation about the share / gj axis is completed. 




r • A ^ 

cos a,, sina^ 0 

-sin a. cos a, 0 


0 


0 1 


(186) 


A 2 rotation is now completed about the gj axis (Figure 58c) to create the 
orientation skew angle jS. 


*2(i3) = 


COsP 

0 

-sinj3 

0 

1 

0 

sin/3 

0 

COsP 
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(187) 







The distinct difference between reaction wheels and CMGs occurs in the third, 
gimbal rotation. This is where CMGs derive most of their advantage. The final rotation is 
a 3 rotation by angle 5 about the new axis (Figure 58d). For reaction wheels, it is 
assumed that 5 = 0. 



cos 5; 
-sin 5; 


V 


0 


sin 5,. 0 

cos 5; 0 


(188) 


The previous series of rotations can now be combined into a single transformation 
converting a vector from the body frame to the gimbal frame. 





s5. 

0 


1 

0 

-sp ■ 


cat 

sa. 

0 

-s5, 

c5, 

0 


0 

1 

0 


-sa, 

ca, 

0 

0 

0 

1 

3 

_ sp 

0 

cP _ 

2 

0 

0 

1 


Gi 


cafficd^ - sa^sSi 
—sa-c5^ — ca^fisd- 
ca^sP 


ca^sd^ + 

caf5. — sa^cfisd- 
sa^sP 


-sficd, 

sfis5. 

cf5 


(189) 


It follows from Equation (45) that is now the transpose of , which is 
equivalent to the reverse sequence via the negative angle of each rotation. 


ca^PcS- - sa-s5. —sa^cS^ — cafPs5. 

-|G 

ca.sP 

ca^s5i + safPc5^ 

cUicS^ — safPsS^ 

sa^sP 

-sped. 

spsd, 

cp 


(190) 


As an example, place the initial CMG at angle a = 0°. The resultant 
transformation from the G frame to the B frame is shown in Equation (191). 
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cf5cd 

s5 

-sj5c5 


-cfisd sP 
c5 0 
sl3s5 cj8 


( 191 ) 


An array of four CMGs is typically arranged in a square pyramid (similar to 
reaction wheels) in whieh eaeh of the CMGs are placed 90° apart. The arrangement used 
here will be; ctj = 0°,a2 =90°,a3 = 180°,a4 = 270°. Figure 59 demonstrates this 
placement. This orientation matches that shown in [29, p. 3] given that the author’s initial 
gimbal angles were each set at 90° rather than d. = 0° as assumed in this thesis. 


CMG4 



CMG 1 


CMG 2 


Figure 59. Body Frame to CMG Sehematic 


Based on the wheel axis deseription mentioned above, the gimbal and wheel 
angular momenta in the gimbal frames ean be written as 
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( 192 ) 


K; = 

1 

o o 

_1 

Gi 

and h^[ = 

1 

0 


1 


0 


Transforming angular momentum from the G frame to the B frame using (184) 
and a = 0° yields Equation (193). Note that the gimbal momentum in the body frame, 
Equation (193), is equal to the reaction wheel momentum in the body frame. Equation 
(53). This is because the orientation of the gimbal rotation axis in is equal to that of the 
reaction wheel spin axis, and both were subject to the same two initial rotations. 



1 

_1 

B 

1 

1_ 

II 

0 

cp 

and = 

1 

1 

1_ 


or in general form: 


- 

B 


- 

k-l. 


and hl^ = 

4 * 1 ,, 

4 * 1,3 

K 

k-k 

4*^,3 



k-k 


(193) 


(194) 


Recall the total angular momentum of the system is the sum of each of the MEDs 
in the body frame. Unlike the single momentum set for reaction wheels in Equation (69), 
CMGs derive their momentum from two different transformation sets. This allows the 
gimbal momentum and gimbal torques to be translated separately to through each of the n 
CMGs to the body frame. 


B 

” 


- 

G 

1 

2® 

_1 


4*kl3 

k'K 

k'L, 


II 

4*,?l3 

k^gi, • 

•• 4*^3 


4. 

i=l 

4*,?l3 

k'gi. 

k^L, 


4„ 
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B 

and = 


where 


and 



k'Si, 


k'Si, 



Br 



k'k 

k • k. 

k'k. 

k'k. 

h'k 

k'k. 

k'k 

k'k, 


k'k, 

k'k 

h'k, 


ys. 

7 ^ 

G 

h 

brSn, 


‘-H-l 



h 

k'k, 


W2 

h'k, 


h 





k'Sn, 

k'L, 

h'L, 

k'k, 

k'k, 

k'L, 


(196) 


(197) 


(198) 


Now, Equations (195) and (196) can be rewritten as 




z,3 


(199) 


K = I.K, 


‘^zfh 


G 

w 




1 

2? 

_1 

G 

1 

1_ 

where 


K 

and = 

K, 



h 


h 



8n 


tv„ 


( 200 ) 


( 201 ) 


For the square pyramid mentioned, the orientation matrices is shown in Equation 
(60) for CMGs. 
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\G 


( 202 ) 



B, MODELING CMG ORIENTATION MATRICES 

These CMG orientation matriees ean be ereated within Simulink. Figure 60 
demonstrates how Equation (190) is implemented. Sinee a, and P are fixed for the 
CMG, their respeetive rotation matriees will be held eonstant. The CMG transformation 
matrix is only state dependent upon 5^. The eolumn seleetors simply extraet the 
first, seeond, and third eolumns individually. 


CMG Orientation Matrix Column 



Figure 60. CMG ith Orientation Matrix Column Model 

Additionally, Figure 61 illustrates how the and R 2 matriees ean be ereated for 
a generie angle. 
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Figure 61. Rotation Matrix Model 


Figure 62 fully ereates all three orientation matriees. from Figure 60 is 

represented four times, onee for each CMG. The individual columns are independently 
combined to create the appropriate 3x4 orientation matrices . 


no 





Orientation Matrices 


B^GI 






Figure 62. CMG Orientation Matrix Model 


C. DIFFERENTIAL EQUATIONS OF MOTION 

The momentum of a CMG ean be defined as the sum of the wheel angular 
momentum and the gimbal angular aeeeleration. 


h,,^=h:+hi 


Each the wheel and gimbal can be independently modeled. 


TJ T N 


T ( N S l B 


/ N B . r tS U, I X 


G, 


(204) 


(205) 
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1. Gimbal Momentum and Torque 

a, CMG Gimbal Momentum - Single Gimbal 

The gimbal momentum is shown first in matrix form. 

hin 0 0 0 

_ 0 0 \ [ 

Shown in dyadic form; 

^gi ~ {jginSiiSii jgi338i38i3)* (^i.?;3) 

~ jghi^iSi^ (‘^'3 * 813} 

~ jgh3^‘Si3 

Change of basis to the body frame; 

K = jgJiSis • + kk+kk) 

= Jgji (^.3 • k )k+jgiji ( 1,3 •k)k + jgji ( 1,3 • k)k 

The resultant rotation matrix is shown as 



(208) 


(209) 


( 210 ) 
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b. CMG Gimbal Momentum - Gimbal Array 

Similar to reaction wheels, the resultant orientation matrix can be created 
using the components from the column of the gimbal frame rotation matrix. 



where 



(213) 


(214) 
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c. CMG Gimbal Torque 

The torque of the gimbal ean now be determined with respeet to the body 
frame. Sinee this is simply the gimbal axis, the orientation is fixed and Z = 0 . The inertia 
is assumed not to ehange relative to the gimbal axis. 




dt 


dt 


(215) 




BrwG 




(216) 


In dyadic form, a single CMG’s gimbal torque can be found using the 
derivative of Equation (210). 


dt 




(217) 


The inertia of the gimbal is not changing with respect to the gimbal frame. 
Also, the orientation of the 3 axis of gimbal frame is not moving with respect to the body 
frame as determined in Equation (183). 

433 =^ andg,^=0. (218) 


Equation (217) now becomes 


dt 



(219) 


( 220 ) 
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Wheel Momentum and Torque 

a, CMG Wheel Momentum - Single Wheel 

The CMG wheel momentum is shown first in matrix form. 


iL G, 7 B G: I 7 G, - - W, 

K, - Jwi + Jw, 


1 

o 

o 

_1 

G, 

0 

G, 

1 

O 

o 

_1 

G, 

’ ’ 

o 

o 


0 

+ 

o 

o 


0 

1 

o 

o 

1_ 


<5, _ 


1 

o 

o 

1_ 


0 


7 Q 

•/ W,’ j j i 

0 


An equivalent derivation is shown in dyadic form. 

K‘ = ) 

= • I, ) + (^,-3 • h ) 

Change of basis to the body frame: 

K ^ ) * (^ 1^1 + kk + kk ) 

= (l,, • )^i+(l,. • k )k +• k)k 

+j.,Ji(gi,*k)k+J.,,Ji(gi, •k)k + pji(h 'k)k 


( 221 ) 


( 222 ) 


(223) 


The resultant rotation matrix is shown as 
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Gi 

p -| 

brki, 

brgi. 

brki. 


j^nPi 

k-k 

k'8i. 

k'h. 


0 

k-l, 

k’l. 

k’8t. 


• 6 

^ tVj 22 ^ 


[k^h)jw,,P-i + [b3^8i)j^,Ji 


( 224 ) 



1 

• 

_1 

B 

k'k 

K = 

k-h 


k’h 


k-k 


1 

• 

Oq > 

1_ 


(225) 


b. CMG Wheel Momentum - Wheel Array 

The resultant orientation matrix can be created using the components from 
the and 3‘^‘* columns of the gimbal frame rotation matrix. 

f r -\B r \ 




k*h, 


k’8i. 


K=t 

: _ 1 


k-k 


k'l. 

jWj 32 

1 — 1 


k-l, 


kPi, 



1 



1- ^ 

) 


k'p 

k'p 

b •£ 

^1 <5rtj 
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1 

J 

P 

_1 

k-p 

k' 82 , • 

•• k'p 


jw2 ip2 

k-p 

4*^2, 

k'L, 


7 Q 

•/ w„ j 1 n 


B 

- 


- 

G 

1 

_1 


Vft3 

k'K 



VVJ 33 i 






j 5: 

+ 

^2*^13 

^2*^23 • 

•• ^2*^„3 


J 1^2 33 1 


4*^3 




• 6 

33 ^ 


(226) 
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iB B fyG f 6G I B rjG f /~\G 

K = ^3 + Zj 


( 227 ) 


c. CMG Wheel Torque - Wheel Array 

The torque of the gimbal ean now be determined with respeet to the body 
frame. Sinee this is simply the gimbal axis, the orientation is fixed and *Z^ = 0. The 
inertia is assumed not to change relative to the gimbal axis. Also, the wheel is kept at a 
nominally constant speed. 


T 


B 

W 



(228) 


(229) 


Using partial differentiation of the second term within Equation (229), the 
following is obtained. 




BryG 


3(5 dt 
dd 3 
dt 3(5 






Therefore, 


where 



can be interpreted as a CMG Jacobian matrix. 


cafj5c5^ - sa^sS- 

-saf5j - ca.cf5s5. 

ca^sP 

ca^sS^ + safficd^ 

caf5j — sUjcPsdj 

sa^sP 

-s[3c5j 

sf5s5j 

cP, 


(230) 


(231) 


( 232 ) 
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By inspection of the transformation matrix, the partial derivative of the 
first column is equal to the second column. 


V ‘ / 2 


(233) 


Similarly, this is captured in the dynamic orientation matrix [6, pp. 464- 
467], where ^ is introduced as a steering logic matrix. 


ft = A5h? where A ■■ 


(234) 


Upon close inspection, A matrix is equivalent to the second column of the 


Z matrix utilized within this thesis. 


A = Zf 


(235) 


An alternative proof follows using dyadics. 


7 ] Q 




(236) 


Using gj = dg 2 from Equation (183), Equation (236) becomes 


h • I, 




(237) 


In summary, all variations are equivalent 

= Aj,a‘s = ^{'’z^)j,n‘s 

"Z” = ‘zfi = AS = —(‘zf)S 

‘ " dS^ ' ^ 


(238) 
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The total wheel torque from Equation (229) now becomes 





( 239 ) 


Assuming like-sized wheels such that Qj = 

becomes a scalar term, the equation can be simplified into 


'“’ll 


BryG CG , • 

+J 




G BrjG 


“’ll ^ 


zr 5 


(240) 


Equation (240) shows the final torque equation for the wheels of a CMG 
array. The equation is dominated by the second term [29]. The first term is dependent 
upon gimbal accelerations, which are generally slow [27]. The second term is dependent 
upon a slow gimbal angular velocity, but it is amplified significantly by the very high 
wheel rate. This is the “torque amplification” that the literature suggests [27, p. 1] 

[28, p. 160]. 

In dyadic form, a single CMG’s wheel torque can be found using the time 
derivative of Equation (223). 


dt 


{K)= +'Li,p-i[hx'k)k 

+ jwi,p-i (k,i •k)k +(ii • k )k +(kii • k)k 
+ (k,i •k)k+ (k,i •k)k+ (k,i •k)k 
+Liji (k,3 • )^i +Liji (k,3 • k )k+Liji (k,3 • k )k 
+ Jmji (k;3 •k)k+J«iji (ki3 • k )k+Jwiji (k,3 • k )k 
+ Jwiji (^3 • k)k + j«iji (k,3 • k)k + j^iji (k,-3 • k)k 


(241) 


The inertia of the wheel is not changing with respect to the gimbal frame. 
Also, the orientation of the 3 axis of gimbal frame is not moving with respect to the body 
frame as determined in Equation (183). The wheel speed is assumed constant. 


= 0 . iwG, = 0 , hi = khi . ^3 = 0, Q,. = 0 


(242) 


Equation (241) now becomes 
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( 243 ) 


‘'d 


* 4)^3 +* 4)4 + •k]k 
+ j.iji (^3 * 4 ) 4 + j^iji (^3 • 4 )4 + j^iji (^,3 • 4 )4 


V 


dt 




1 

• 

_ 1 

B 

1 

J^> 

• 

_1 

k'Si, 

i ad.+ 

dWj 23 11 

4*4 

4-4 


1 

u> 

• 

1_ 


4 33 4 


The first term again validates the claim made in Equation (237). 

3, Total Momentum and Torque 

The combined CMG gimbal and wheel momentum and torque are now 


KMa=K + K 


=Czf 7 Al+("ztM” + "z»/..£!°) 

= 'z?( 4 ,+ 2 .,,)«“ + 'z,%fl,“ 


^CMG ~^s 


V ^•>833 g I V 3 ./W33''g 2 ./w,, g''J 

= ''z^(j +j )s^ + ‘'z^(j a^)d 

3 \'^8i3 •^'^33 / .? 2 g J 

= {j +j )'^z^d^ + {j n^)’^z^d 

\Jg 33 Jw^-i j 3 y./w,, g j 2^ 

By Euler’s transport theorem, the inertial derivative is 




dt 


dt 

i.B , N 


CMG 


"'CMG ~ "CMG ^ W Hf^MC 


(244) 


(245) 


(246) 


(247) 


Equation (247) can be rewritten by substituting the angular momentum terms 

= (/ + / Yz’^S'^ + ij Q.^Vz'dS 

8 833 ^^33 j ^8 8 ] 2 


+ X 


ffj +;• Yz^5^ + j ’^Z^d^]' 

\V‘'^33 •^^33) 3 g Jw;, 3 g I 
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(248) 








Like reaction wheels, CMGs are subject to the same pseudo-inverse law. Much 
like the pseudo-inverse of the reaction wheel orientation matrix, (or A) can be 

inverted using the same method. Here will be used to map the three body axes back to 
n CMGs. 


= (/ +j )^z^S^ + (j 

req \J gii ^ S «/ ^ 


+ X 




B 


^ 30 ,) 


(249) 


Although Equation (249) shows the full dynamics, the gyroscopic coupling of the 
body motion with the gimbal rate can be assumed as very small [29]. Additionally the 
gimbal acceleration can be assumed small [29]. Therefore, Equation (249) can be reduced 
to Equation (250) and rearranged per Equation (251) using the pseudo-inverse. 

Q.^Yz^5 

'‘req g J 2'' 

= (250) 

BrjG^G 

— Zo T 

2 w 

=(ft°r[”Zt ]*(-<) <251) 


This is the same steering logic proposed in [6, Algorithm 7.147]. Therefore, the 
assumptions above were taken as well. 


D, MODELING CMGS 

In order to apply the physical actuator limitations of the CMG as part of the 
simulation model, the torque required of the assembly must be broken down, individual 
gimbal and wheel limit saturations applied, and then the torque of the CMG assembly 
rebuilt. Within Eigure 63, (251) is used to break down the torque required and allocate it 
across the four wheels. The orientation matrix Z block consists of Eigures 60-62. 
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CMG Steering Law 

(saturation occurs in red) 



This CMG model can be further expanded to include null motion in the presence 
of near-singular conditions. Additionally, the components of Equation (249) can be 
modeled, but Equation (251) is modeled here for simplicity. 

E. CMG MOMENTUM SPACE 

The boundary created by the maximum momentum in any given direction for a 
particular CMG configuration is much more difficult to define than that of a reaction 
wheel configuration. The momentum space is no longer a three-dimensional geometric 
envelope with all interior points possible. The physical limitation is still the vector sum of 
the maximum capacitates of the individual CMGs, where momentum saturation occurs. 
Unlike reaction wheels, the surface is not a polyhedron. Instead, the CMG momentum 
envelope is contains singular cavities that pass through the center from one edge of the 
momentum space to the other. 

Eike the reaction wheel case, maximum momentum space for an evenly spaced 
(Act = 90) four-CMG array is when = 54.73° [30, p. 866 ]. The resulting surface is 
shown in Eigure 64. Since the gimbal axis ^3 and the wheel angular momentum g, are 

always perpendicular, a reduced amount of wheel momentum is available on the surface 

122 


































of the momentum envelope, which is set within the body frame 

= ^ hr^cn,gA = Si,’kmgA = '^]- Thcsc Tcgions are identified by white 

cavities in Figure 64, in which 2" cavities exist [28, p. 161]. The identification of 
singular surfaces will be discussed in Section E. 



Figure 64. Momentum Space of a 4 CMG System Given p=54.73° 


The shape within Figure 64 is unique to this four-CMG system. The shape alone 
reveals configurations of maximum momentum. The momentum peaks help identify the 
best possible axes to rotate about. These maximums occur at = 3.266/z^^, revealing 
that when using CMGs in comparison with similarly sized reaction wheels, 40% more 
momentum is achievable about certain axes. 


123 










The concept of placing a momentum sphere directly inside the minimum radius of 
this shape does not make sense because the shape is concave. This leads the discussion to 
CMG singularities. 


F. CMG SINGULARITIES 

Within a singular gimbal angle configuration, it is not possible to produce a 
torque in certain directions. Therefore, the required torque may or may not be produced. 
Singular states occur when the rank of is not full (less than n CMGs). This condition 
can occur because is state dependent. 

B r -iG, 

-sa.c5. - ca.cPsS. 

(5;) = ca^cd^ - safPsd^ (253) 

sps5, 

If a configuration of delta leads two CMGs to yield the same columns within their 
respective Z matrix, then a singular state occurs. The entire array for the [0,90,180,270] 
arrangement is shown in Equation (254). 

Br -iG 





-cd^ 


c5^ 



c(5s5^ 

(254) 

s^sd^ 

s^sd^ 


sf5s5^ 



For example, when = 0° and = 180°, two columns are the same. Therefore, 
rank(^Z^) = 2, which is not full. Therefore, a singular condition results during the 
pseudo-inverse. 

B - 

0 ^^^2 0 

1 -cpsd^ 1 

0 s(5s52 0 
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This is just one of many different eonfigurations ean reduee the rank and eause a 
singular eondition. 


1. Solving Singularity Conditions 

A detailed proeedure and set of equations for determining the singularity 
conditions and associated singular momentum vectors is given in [6, pp. 676-678]. CMG 
singularity occurs when the magnitude of the dot product of the wheel momentum vector 
and the spacecraft rotation axis is one: g. • = hf • = ±1 • In this thesis, only the 

results of the four CMG case will be analyzed. The singularity condition for each CMG 
can be written as 

h ^ ^cmsA ) X K _ ^ (g-3 ^ ^ k 

gi.Xh^mgA 




(256) 


e, 


h^=±- 


The vectors e, contain the cross product of the individual CMG gimbal vectors 

1 B 

and the total body momentum unit vector h^^gA ■ Both positive and negative variations are 
accounted for within Equation (257) by the sign function of Equation (258). 


1 - ( + Cp kmgA,^ ) 

G = ± jl - ( kmgA,. + CP kmgA, ) 


(257) 


G = ± - (-^/5 kmgA, + Cp kmgA^ ) 


e, = +J1- - 


■ -SpkmgA,+cfikmgA) 


e^=sign{e.) = ±\ (258) 

The total CMG array momentum in the body frame is created using Equation 


( 259 ). 


KmgA=J,hf = J^ 




(259) 


«=1 


«=1 
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Individually accounting for each body axis separately, Equation (259) can be 
written as Equation (260). 




bf = 


, c 


B 

cmgA^ 




h,= 


P(sPh^n,gAr^P^cmgA^) ^ ) 


(260) 


+ ■ 


P ( sP h^gA + cP + cp ) 


Given four CMGs, there are 16 (2") combinations of q that produce different 

singular conditions. When q has the same sign for all CMGs |e = ±[l 1 1 l]j’ 

the outer singularity surfaee shown in Eigure 64 is produced. Other eombinations result in 
interior singularities: states that are eontained within the volume of the CMG momentum 
envelope. 

Interior singularities mean that at some non-maximum momentum (within the set 
defined by the geometry of Eigure 64), the required torque magnitude or direction cannot 
be ereated exactly, creating an off-eigenaxis rotation. Two types of internal singularities 
exist: hyperbolic and elliptic singularities. The names derive from the shape of the gimbal 
angle motion in the neighborhood of a singular point. These are discussed in detail within 
[3, pp. 23-30]. 


2. Hyperbolic Singularities 

Hyperbolic singularities are the milder of the two internal singularities, 
considered passible in [3, pp. 23-30]. This is the case when the CMG array cannot 
exactly produce the required torque. A speeific gimbal signal ean be added to the steering 
logie whenever a singular state is detected. This null motion produces no net torque on 
the system, but assists in passing through hyperbolic singularities. They oceur when 
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exactly two of the e. have mismatched signs. The six combinations are shown in the 
following three cases. 


e ,=± 

1 

1 

-1 

> e,- = ± 

1 

-1 

1 

’ £/ = ± 

1 

-1 

-1 


-1 


-1 


1 


(261) 


The hyperbolic singularities make up the very center of the shape, extending from 
the origin (zero momentum) to a maximum momentum magnitude of 1.897 at its largest. 

A graphical depiction of case 2 |e = ±[ 1 -1 1 -1 ] | is shown in Figure 65. It is 

noteworthy that it is symmetric about all three body axes. 



•05 


Hyperbolic 
Singularity 
(case 2) 


f. = ± 


16 










Figure 65. Hyperbolic Singularity Surface for 8i=±[l -1 1 -1]^ 
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Figure 66 demonstrates all three cases plotted together. Each of the cases 
extended from the origin to 1.897 . 



Figure 66. Flyperbolic Singularities 

3, Elliptic Singularities 

Elliptic Singularities occur when the torques allowable by the physical layout of 
the CMGs are perpendicular to the required torque [11, p. 179]. They are impassible [3]. 
Elliptic singularities occur when only one e. sign is mismatched. The eight possible 
combinations are shown. 
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The elliptic singularities extend across the width of the shape, extending from a 
minimum of 1.0069 to a maximum of 2.99Sh^^^ at its largest. A graphical depiction of 

easel |e = ±|^-l 1 1 l]| is shown in Figure 67. It is noteworthy that the singular 

surface is oriented at the skew angle (^ = 54.73°) due to the alignment of the CMG. 



Figure 67. Elliptic Singularity Surface (case 1) 


Given all cases at once, independent plots extend across the momentum space, 
shown in Figure 68. 
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Figure 68. Elliptic Singularities 
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Most interestingly, the eight interior horns shown in Figure 68 align with the eight 
hollow circles shown in the perimeter volume in Figure 64. This union is shown in Figure 
69. 
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Figure 69. Elliptic and Surface Singularities 


G. CMG LARGEST INSCRIBED SPHERE 

This thesis will not focus on complicated CMG singularity avoidance algorithms, 
an ongoing and developing research area in its own regard. In CMG systems, it may be 
easiest to treat these singularities simply as impassible in order to avoid exotic 
controllers. To create a volume that is guaranteed reachable, the minimum momentum 
magnitude of the elliptic singularity region must be completely avoided. Since this 
minimum occurs at the calculated values of 1.0069/z^^^ (1-0 ^max this also creates 

the largest sphere of guaranteed momentum capability. Superimposing a sphere into the 
elliptic region creates the space of guaranteed mobility, shown in Figure 72. 
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Largest Guaranteed Passible Sphere 



Figure 70. CMG Largest Sphere 


This reveals that the maximum guaranteed momentum capacity is 1.0069 , 

69% smaller than originally estimated in Equation (252). 


hi 


h. 


1.0069h„ 

3266h. 


= 0.3083 


(263) 


This means the available momentum in the CMG case is actually less than that of 
the reaction wheel cases: recall 1.633 for the full reaction wheel and 1.333 for 

IIlaA IIlaA 

the pseudo-inverse limited reaction wheel. 


H. SUMMARY 

This chapter was used to develop the CMG equations of motion from the first 

principle of conservation of angular momentum. The additional moving frame causes 

CMG dynamics to be much more complex than reaction wheel dynamics. The orientation 

matrix was shown to be state-dependent, which leads to singularities. Certain CMG 

arrangements lead to different types of singularities that can cause controller failure. 

Elliptic singularities are non-passable and hyperbolic singularities are passable. 
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The singularities were individually mapped in three-dimensional momentum 
spaee to demonstrate where they oeeur. The momentum envelope was eombined with 
both of the interior hyperbolie and elliptie singularities to demonstrate their interseetion. 
Using the simplifieation that elliptie singularities must be avoided and hyperbolie 
singularities are passable, a sphere of aeeessible momentum was identified. The radius of 
this sphere guarantees the momentum of 1.0069/Zj^^ about any axis of rotation. 
Therefore, this momentum value ean be used to maneuver the CMG about any eigenaxis. 

Despite the lost momentum capaeity and required eomplex ACS architeeture, 
torque amplifieation remains the primary driver for CMG desirability over reaetion 
wheels. For CMGs, a small gimbal input results in a large torque output. This relationship 
is quantified in Chapter VII. 
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VII. REACTION WHEEL AND CMC COMPARISON 


The previous ehapters were used to develop the neeessary baekground and 
equations of motion to support a detailed diseussion within this ehapter. This ehapter is 
used to help bridge the performance gap between reaction wheels and CMGs. Instead of 
insisting that CMGs must be used for highly agile spacecraft, this chapter will introduce 
some concepts for sizing reaction wheels to attain the required slew performance. In 
some cases CMGs will be more desirable, but this chapter will illustrate that for some 
satellite systems, reaction wheels can potentially outperform CMGs. 


A. MOMENTUM SPACE AS A VISUALIZATION TOOL 


In Chapter IV, the reaction wheel momentum space was discussed. Due to the use 
of the pseudo-inverse for control allocation, the entire momentum and torque space was 
not available for control. Also, there is no real reason to stay on the eigenaxis, despite the 
fact the eigenaxis maneuvers have been historically implemented on spacecraft. The input 
can be shaped in quaternion space to move around the momentum space, searching for 
the axes of maximum momentum and torque (vertices from Chapter IV). This allows the 
slew time to be reduced and/or optimized, regardless of the slew distance and path 
traveled. Recall from Chapter V that for a successful maneuver, it is the boundary 
conditions that must be satisfied, not the path itself 


Consider a 120° eigenaxis slew about e = 


yfin 0 2/V3 


T 

. The same 


spacecraft/reaction wheel configuration from Chapter V, Sections F-G will be used. 
A trapezoidal input is fed into the system, accounting for the limitations of the pseudo¬ 
inverse. The dynamic response of the reaction wheel array is used to build Figures 71-72. 
The vectors in the body frame are normalized against the individual wheel maximum. 




I* 

norm 


T It 

•' ' • ma\ 


(264) 


For a reaction wheel system, the shapes of the momentum and torque envelopes 


are exactly the same. Thus, normalization via (264) allows both momentum and torque 
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signals to be compared on the same scale. Figure 73 shows the scaled system response to 
the shaped input. Only the relative magnitude is shown with no regard for direction. The 
system follows the signal exactly: bang-coast-bang. The only deviation is a discontinuity 
at Point 5. This occurred because the shaped trajectory implemented cut the velocity off 
when the angle reached 0.999 of the final angle, causing a split second of unpredictable 
acceleration, and is an artifact of the simulated model only. 



Momentum 

Torque 


Time (sec) 


Figure 71. Normalized Momentum and Torque for a Trapezoidal Response 


The same maneuver is next mapped onto the momentum and torque space in Figure 
72. Since the maneuver is an eigenaxis maneuver, the momentum travels from zero to 
maximum momentum within the envelope defined by the pseudo-inverse control 
allocation, yet the axis of rotation remains fixed. The numbered points reflect those from 
Figure 71. On the right is the torque mapping. The locations of Points 1,2 are maximum 
torque in the direction of the eigenaxis. Point 5 demonstrates a strong torque reversal to 
decelerate the system. 
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regions mathematically unavailable 


Figure 72. Momentum/Torque Space During Eigenaxis Slew 


The results shown in Figure 72 are uninteresting, and it is not surprising that this 
envelope mapping is not used to help visualize conventional reorientations. The same 
maneuver will now be illustrated using the optimal control approach described in Chapter 
V, Figures 51-53. 
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Figure 73. Normalized Momentum and Torque for an Optimal Maneuver 


Now the temporal element of the identified points becomes important. The following 

is a description of the instantaneous rotation axis as it pertains to Points 1-5 given in 

Figures 73 and 74. 

Point 1. The array pulled far off eigenaxis at the point of maximum torque, a maximum 
vertex at 2.309 . This allowed the system to gain momentum as quickly as 

possible. 

Point 2. The array quickly moved along the edge of the torque space to a minor vertex at 

2.000 T_. 

Point 3. The array rapidly changed directions again, this time to another maximum vertex 
at 2.309 . Next, the torque decelerated the system as it pulled the rotation 

angle close to the eigenaxis. 

Point 4. The rotation angle wobbled around the eigenaxis, forming a lasso-like loop near 
maximum momentum. The torque remained close to zero, making a few 
adjustments necessary to keep the orientation direction on track. 

Point 5. Point 1 is reversed to decelerate the array. 
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Figure 74. Momentum/Torque Space During Optimal Maneuver 


The type of action shown in Figure 74 is common for most optimal maneuvers. The 
activity appears unintuitive, but trends develop. Torque seeks high ground to raise the 
momentum quickly. Torque and momentum follow the edges and vertices as they move 
around. Rarely does the maneuver stay very close to the actual eigenaxis. Figures 75 and 
76 represent two more maneuvers to further illustrate some of the interesting behavior, 
obtained as a result of optimizing the slews. The same momentum limits are used, but the 
torque is doubled to 0.2 Nm to illustrate different maneuvers. 

The optimal use of the full momentum and torque envelopes is extremely interesting. 
The momentum and torque move like water encased in a polyhedron, flowing to high 
surfaces. The system always seems to tend towards maximum torque in a particular 
direction to maximize the momentum in order to rapidly maneuver the spacecraft. 
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Figure 75. Optimal Slew Example 2, qo=[0.36 0 0.79 0.5], qf=[0 0 0 1] 
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Figure 76. Optimal Slew Example 3, qo=[0.23 0.17 0.82 0.5], qf=[0 0 0 1] 
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B, MOMENTUM AND TORQUE RELATIONSHIP 


Many variables go into determining the torque output of reaction wheel and CMG 
systems. This section uses market data to develop a “rule of thumb” relationship that can 
be utilized in size comparisons between the two very different systems. 


1. Simplified Dynamics 

In order to size a MED array, a reasonable approach would be to use the largest 
inertia. In the same way that the minimum momentum limits angular rate, the largest 
inertia limits rate. A starting point for sizing wheels might begin without the 
complications of gyroscopic coupling per Equation (265). 


h 


N 

medA 


T 


N 

medA 



(265) 


In Chapter III, it was shown that the spacecraft can be modeled as a simple double 
integrator. This creates a linear relationship between spacecraft rate/acceleration and 
momentum/torque, as both use the same proportionality factor J. This allows a very basic 
relationship between the dynamics of the spacecraft and the momentum exchange device 
to be considered. As the engineer narrows down selection within a few reaction wheels, 
subsequent iterations should include gyroscopic coupling in the sizing analysis. 


2, Market Research and Simplified Momentum and Torque Equations 

Market data may be used to develop a relationship between momentum and 
torque for reaction wheel and CMG systems. Appendix B consists of a table of 
performance data collected from company specifications sheets and published papers 
across the Internet. Reaction wheels are more common and over fifty data points were 
readily available. CMGs are far less prevalent as less than twenty systems could be 
identified. Figure 77 offers a glimpse at this momentum and torque comparison for 
reaction wheels and CMGs. Reference [4, Figure 4] offers a similar chart and the six 
CMG data points presented in this paper are overlaid on the chart for comparison with the 
information presented in this thesis. The market data from Appendix B are marked with 
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“X’s.” The CMG data from [4] are marked with circles. Common points contain both 
marks, but are counted only once to facilitate curve-fitting of the data. 



0.0001 0.001 0.01 0.1 1 10 100 1000 10000 


Specified Angular Momentum (N-m-s) 

Figure 77. Market Research Torque vs. Angular Momentum for Reaction Wheel and 

CMG Systems 

The data in Figure 77 is curve fit to obtain the trend for each type of system. A 
few things become clear. First, there is a near proportional balance between CMG 
momentum and torque. Equation (266). This is because the individual CMG torque is 
determined by the product of gimbal rate (commonly 1 rad/s [27, p. 8]) and momentum. 
On the other hand, reaction wheels seem to have an approximately square root 
relationship between momentum and torque with a proportionality constant of less than 
three percent. Equation (267). The market data confirms that reaction wheels have far 
less torque capability than CMGs. 

CMG: since S^lrad/s (266) 
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(267) 


Reaction Wheel: ^ 0.027yfhZ 

where torque units are N-m and momentum units are N-m-s. 

Figure 77 is created on a per-wheel/CMG basis. This means that MED arrays are 
not accounted for, so array calculations using orientation matrices are still necessary. 
Another notable figure is the large torque difference at high momentum levels. This 
reveals why so many CMGs are being developed in the 100+ Nms range. On the other 
hand, micro-CMGs are rare, seemingly for good reason. With the shape of the curves, 
there actually may be a torque advantage for reaction wheels around the 0.001 Nms 
momentum level. Regardless, these curves offer a basis for comparison. Therefore, given 
momentum storage of a reaction wheel or CMG, the torque available will be described as 


CMG: 

„ / „ \1.0019 

=1.0967 

Cmg v.' » j 

(268) 

Reaction Wheel: 

w / w \0.5351 

tZ=00270{hZ) 

(269) 


C. COMPARISON DEVELOPMENT 


For comparison purposes, a propagation of effective eigenaxis trajectories for the 
optimal models will be completed, similar to Figures 51 to 53 from Chapter V. The 
optimal maneuvers will be propagated in a bang-bang or bang-coast-bang trajectory 
about its effective eigenaxis (Chapter V, Section H). Already derived previously in this 
thesis, the following spherically-available momenta will be used. 


w 


hi, =1.333h 


max 

w 


- 577 / 1 , 


= 1 - 805 ^^ 


Ca =l-007h^ 


(270) 


The torque equations within Equations (268) and (269) will be used to compute 
reaction wheel and CMG output torques. Transformation of torque available to the body 
frame, given a 4 reaction wheel or 4 CMG array, is computed as follows. 
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.B 

rwA 

= 1.333tI 

.B 

rw-pseoptA 

= 1.333t: 

.B 

rw-optA 

= 1.333t;: 

.B 

cmgA 

= l-007< 

' '^’^max 


(271) 


For 1 Nms wheels, the data in Table 3 ean be determined using Equations (270) 
and (271). 



Wheel/Gimbal Frame 
hmax (Nms) T^ax (Nm) 

Body Frame 

hmax (Nms) Tmax (Nm) 

RW 

1 

0.027 

1.333 

0.036 

RW pseudo-opt 

1 

0.027 

1.577 

0.036 

RWopt 

1 

0.027 

1.807 

0.036 

CMG 

1 

1.097 

1.007 

1.104 


Table 3. Body Momentum/Torque Comparison Table - 1 Nms Momentum 


Now a spaeeeraft size elass must be established. Sinee the MED momentum 
storage is the starting point, either the spaeeeraft size or desired slew rate must be ehosen. 
Spaeeeraft rate will be ehosen next, beeause we are focusing on agile spacecraft. In the 
space industry, a 2-3 °/s slew is considered agile, even for CMGs. Eor example. 
Worldview 1 slews at 2.3°/s [2]. Therefore, a rate of = 3° I s will be the benchmark 
with which the spacecraft sizing will be determined. The base reaction wheel will be used 
for sizing, and all other systems will be required to slew the same spacecraft. 


1 B 1 B ■ N ^^B 

^rwA ~ ^req ~ J sic ® 


rwA 


‘ J sic 


N ..B 

(O 


1.333 Nms 
3° Is 


= 25.5 kg-m^ 


(272) 


Using these properties, the individual spacecraft rates for the various systems can 
be solved for. 
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(273) 


(274) 


The data from equations (268) through (274) are eompiled in Table 4. Notiee the 
higher momentum eapacity for reaetion wheels eorresponds direetly to the larger 
maximum spaeeeraft rotation rates. Note that this rotation rate only refers to the effeetive 
rotation rate about the eigenaxis. The aetual rotation rate of the optimal ease about the 
instantaneous axis of rotation ean potentially be mueh faster so the values in Table 4 are 
eonservative estimates. The body aeeeleration of the reaetion wheels is more than an 
order of magnitude lower than the aeeeleration possible using the CMGs. 



Wheel/Gimbal Frame 

hmax(Nms) Tmax(mNm) 

S/C 

js/c (kg-m^) 

Body Frame 

hmax (Nms) Tmax (Nm) 0 ) (7s) a (7s^) 

RW 

1 0.027 

25.458 

1.333 0.036 3.000 0.081 

RW 

’ ’ pseudo-opt 

1 0.027 

25.458 

1.577 0.036 3.549 * 0.081 * 

RW„p, 

1 0.027 

25.458 

1.807 0.036 4.067 * 0.081 * 

CMG 

1 1.097 

25.458 

1.007 1.104 2.266 2.485 

*effective rotation and acceleration rates given about an effective eigenaxis. Refer back to Section H of Chapter 6 for 
description of effective eigenaxis. The actual body rates and acceleration can be much higher throughout the maneuver. 


Table 4. Body Rate/Aeceleration Comparison for a 25.5 kg-m^ Class Spaeeeraft 
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Given the effective body rate and acceleration data from Table 4, the data can 
now be propagated through the “effective” or equivalent eigenaxis slews using 
trapezoidal input trajectories. This is the same as what was done for the optimal reaction 
wheel case from Chapter V. Recall the curve from Chapter V was built using a specific 
reaction wheel of 0.1 Nm torque, much higher than the one used here per the relationship 
in Equation (269). Therefore, the curves are built using different acceleration data, but the 
concept is the same. 

The plot of the eigenaxis slew angle versus slew time for all four cases is shown 
in Figure 78. This figure demonstrates the shape of the curves. Since the reaction wheel 
arrays have higher momentum storage, the time/angle slope of the momentum-limited 
region is much smaller than the one for the CMG case. This causes the reaction wheel 
systems to eventually catch up with the CMG system. The time-optimal case would not 
catch up with the CMG until the maneuver exceeds 180°, which is an unrealistic slew. 
For a 25.5 kg-m^ class spacecraft, the control system comparison shows CMGs to be 
more favorable for agile maneuvering. 
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Now, suppose the wheel momentum is deereased by an order of magnitude 
(^max =0.1A^ms). Using the same = 3° I s requirement, the corresponding spacecraft 
is now 2.55 kg-m^. The new comparison data is shown in Table 5. 



Wheel/Gimbal Frame 

hmax(Nms) Tmax(Nm) 

S/C 

js/c (kg-m^) 

Body Frame 

hmax (Nms) Tmax (Nm) 0 ) (7s) 

a (7s") 

RW 

0.1 

0.008 

2.546 

0.133 0.010 

3.000 

0.236 

RW 

’ ’ pseudo-opt 

0.1 

0.008 

2.546 

0.158 0.010 

3.549 * 

0.236 * 

RW„p, 

0.1 

0.008 

2.546 

0.181 0.010 

4.067 * 

0.236 * 

CMG 

0.1 

0.109 

2.546 

0.101 0.110 

2.266 

2.474 

*effective rotation and acceleration rates given about an effective eigenaxis. Refer back to Section H of Chapter 6 for 
description of effective eigenaxis. The actual body rates and acceleration can be much higher throughout the maneuver. 


Table 5. Body Rate/Acceleration Comparison for a 2.55 kg-m^ Class Spacecraft 


Notice the spacecraft rates did not change across the two models. Instead, the 
spacecraft acceleration is the important change. Due to the roughly 1:1 ratio CMC 
momentum to torque ratio ~ from Equation (268), the CMC trajectory did not 

change very much. On the other hand, the reaction wheel has a square root relationship, 
~ 0.027^/^ ) from (269). This causes the reaction wheel torque to momentum ratio 

to improve as the wheel size is decreased. The reaction wheel acceleration rates are now 
about three times higher than for the larger wheel (see Table 4). Generating the effective 
eigenaxis trajectories for each system gives Figure 79. Using these curves, it becomes 
evident that any slews greater than 77° will favor the time-optimal reaction wheel control 
system over the CMG. 
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Notice the significant change in the reaction wheel curves. The spacecraft rates 
did not change across the two models, but the improved/larger spacecraft acceleration 
obtained from the reaction wheel system increased by about a factor of three. Because the 
angular rates are the same as before, the time/angle slope of the momentum-limited 
region did not change at all. Due to the higher available torque, the momentum-limited 
region starts at a lower angle for the reaction wheels. This causes an earlier intersection 
with the CMG curve. To better illustrate these transition points, the reaction wheel curves 
are looked at as a function of angle, after which the performance of the reaction wheel 
assembly exceeds that of the CMG assembly. This threshold will be referred to for the 
remainder of the thesis as the reaction wheel advantage angle, shown in Figure 80. 
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Figure 80. Maneuver Angles where Reaction Wheel Systems Outperform CMGs for a 

2.55 kg-m^ Class Spacecraft 

The trend of rapidly increasing reaction wheel performance is continued when 
wheel momentum is decreased by another order of magnitude {h^^=0.0lNms). 

Using the same ^(0^ = 3° Is requirement, the corresponding spacecraft size is now 
0.255 kg-m^. The comparison is shown in Table 6 and Figures 83-84, illustrating that 
the reaction wheel advantage angles decrease as body accelerations provided by the 
wheels increase. 



Wheel/Gimbal Frame 
(Nms) Tmax(mNm) 

S/C 

js/c ( kg - m ^) 

Body Frame 

h„a« (Nms) T^ax (Nm) w (“/s) 

Ct ( 7 s ") 

RW 

0.01 

0.002 

0.255 

0.013 0.003 

3.000 

0.689 

RW 

’ ’ pseudo-opt 

0.01 

0.002 

0.255 

0.016 0.003 

3.549 * 

0.689 * 

RW„,. 

0.01 

0.002 

0.255 

0.018 0.003 

4.067 * 

0.689 * 

CMC 

0.01 

0.011 

0.255 

0.010 0.011 

2.266 

2.464 

*effective rotation and acceleration rates given about an effective eigenaxis. Refer back to Section H of Chapter 6 for 
description of effective eigenaxis. The actual body rates and acceleration can be much higher throughout the maneuver. 


Table 6. Body Rate/Acceleration Comparison for a 0.255 kg-m^ Class Spacecraft 
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Figure 81. Slew Angle vs. Time Comparison for a 0.255 kg-m^ Class Spacecraft 
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Figure 82. Maneuver Angles where Reaction Wheel Systems Outperform CMGs for a 

0.255 kg-m^ Class Spacecraft 


Referring to Figures 78, 79, and 81, it is apparent that when using like-sized 

reaction wheels and CMGs, the larger momentum capacity of reaction wheels is much 
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more advantageous for small spacecraft than for larger spacecraft. In the case shown in 
Figure 83, the optimal reaction wheel control system becomes more agile than the CMG 
control system for all maneuvers greater than 24.1°. Standard operation of the reaction 
wheel array does not begin to outperform CMGs until the maneuver size exceeds 30.3°. 
If only the reaction wheel advantage angle is plotted against the wheel momenta, the 
curves in Figure 83 are obtained. The curves reveal that reaction wheels overcome CMGs 
at relatively small angles for small rotor momenta. As the spacecraft momentum 
increases, the reaction wheel advantage angle increases rapidly and reaction wheel 
control systems promptly lose their edge. 
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Figure 83. Reaction Wheel Advantage Angle vs. Wheel/Rotor Momentum; Where RWs 

Overcome CMGs (m=3°/s) 


Likewise, if the spacecraft inertia versus slew angle is plotted (assuming ro=3°/s), 
the trend is shown again (Figure 84). For small spacecraft, reaction wheels seem to 
outperform CMGs, especially for large angle maneuvers. 
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Figure 84. Reaction Wheel Advantage Angle vs. Spacecraft Inertia: Where RWs 

Overcome CMGs (m=3°/s) 

Thus far, the results of this chapter show that reaction wheels outperform CMGs 
for small spacecraft. In addition, the two optimal control approaches described 
significantly outperformed the heritage eigenaxis maneuver. For small spacecraft, the 
reaction wheel advantage angle was reduced by only a few degrees by using optimal 
reaction wheel methods (Figure 85). At larger slew angles, 40-60° improvement can be 
achieved by optimal systems over the heritage eigenaxis maneuver (Figure 85). This 
means for large angle slews, the optimal reaction wheel system outperforms the CMG 
system significantly more quickly than if heritage control logic were used. 
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Figure 85. Zoom in; Reaction Wheel Advantage Angle vs. Spacecraft Inertia; Where 

RWs Overcome CMGs (m=3°/s) 

Figure 86 shows these reductions in reaction wheel advantage angles as a 
percentage improvement. The pseudo-inverse optimal method reduces the reaction wheel 
advantage angle by roughly 15% over the heritage eigenaxis shaped trajectory, while the 
full time-optimal solution reduces the reaction wheel advantage angle by 18% across the 
entire spacecraft inertia regime. 
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Figure 86. Optimal Improvement over normal Reaetion Wheel Systems 


D. MOMEMTUM EXCHANGE DEVICE PERFORMANCE AS A FUNCTION 

OF MASS 

In the analysis of Seetion C, equally sized reaetion wheels and CMGs were used, 
i.e., both had the same momentum storage eapaeity. This does not have to be the case. 
Suppose the systems engineer gives a mass restriction for the entire MED array to the 
ADCS engineer. The ADCS engineer can quickly determine whether reaction wheels or 
CMGs would be better given this limitation. In this section, it will be shown how much 
more momentum can come of a reaction wheel array in comparison with CMGs given a 
mass constraint. Figure 87 shows the market research for individual MED mass and the 
output momentum the actuator can achieve. The data was curve fit to develop a general 
trend. 
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Figure 87. Market Research Angular Momentum vs. Mass 


The mechanical complexity of CMGs requires additional mass and volume within 
the spacecraft. Using this knowledge and the relationships from Section C, the full array 
mass and body momentum can be created using the trend lines (Figure 88). This analysis 
shows that for the same array mass, an order of magnitude more momentum storage is 
typically available from reaction wheel arrays than from CMG arrays. This implies that if 
mass is a critical design constraint for the spacecraft, then engineers should take an even 
harder look at reaction wheels. 
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Figure 88. Array Momentum vs. Mass 


Torque output, on the other hand, is better for large mass CMGs. Figure 89 
demonstrates the transition where a CMG array will begin to create more torque than a 
reaction wheel array. This is an important point because it helps illustrate why few CMGs 
have been developed on the small scale. From the market research shown in Figure 87, 
only one CMG was actually smaller than 0.898 kg. 
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Figure 89. Array Torque vs. Mass 


E. SUMMARY 

In this ehapter, it was shown that reaction wheels seem to perform better than 
CMGs for small spacecraft. As the spacecraft size decreases, the reaction wheel 
advantage increases. This increase can be further improved by the use of optimal control 
algorithms. 

The CMG torque advantage is evident for large moments of inertia. On the other 
hand, the torque gap lessens as the spacecraft size (therefore the MED array momentum 
capacity) decreases. Although the gap does not completely close, the ratio of CMG 
torque to reaction wheel torque decreases significantly with spacecraft size. 

Given arrangements of reaction wheels and CMGs with similar individual MED 
momentum, the reaction wheel array will have a larger momentum storage capacity. This 
causes reaction wheel systems to catch up with CMG systems over large angle slews. The 
reaction wheel advantage angle is the slew size in which a reaction wheels become a 
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better alternative than CMGs. For large spacecraft, the reaction wheel advantage angle 
angle is very high, possibly over 180 degrees. 

Time-optimal shaped maneuvers of reaction wheel systems can be used to lower 
this reaction wheel advantage angle by 13 to 21%, depending upon algorithms used and 
spacecraft size. By shaping only the quaternion and body rate trajectories, the 
mathematical limitations of the pseudo-inverse control allocation must still be ascribed 
to. This allows a reduction of the time-optimal reaction wheel advantage angle by 
approximately 15% over the advantage angle of the shaped eigenaxis maneuver. The 
addition of shaped time-optimal wheel torques reduces the advantage angle by 
approximately 18%. 

Considering the mass of the reaction wheels and CMGs, the reaction wheels 
produce a much higher momentum limit than CMGs for a given mass. Interestingly, the 
additional order of magnitude of momentum equates to additional torque in the similar- 
weight MEDs. For MEDs smaller than 0.898 kg, the reaction wheel array will actually 
produce more torque than a CMG array. This means that given a small spacecraft with 
weight restrictions, reaction wheels are the better alternative for both torque and 
momentum. 


158 



VIII. CONCLUSIONS AND FUTURE WORK 


A. CONCLUSIONS 

This thesis is an attempt to bridge the gap between reaction wheels and CMGs, 
particularly for small satellites. Using the software upgrades and relative sizing concepts 
proposed, reaction wheel arrays have the ability to outperform CMGs for some classes of 
spacecraft. The relatively low cost of reaction wheels can be attractive once again. 

Reaction wheels are typically used inefficiently within the spacecraft industry. 
The standard maneuver, a shaped bang-coast-bang maneuver about the eigenaxis, should 
be reconsidered in the presence of new control algorithms. The results of this thesis show 
that the speed of this eigenaxis-restricted maneuver can be improved 18% in each of 
momentum and torque by applying the pseudo-inverse limited time-optimal shaped 
quaternion and rate trajectories. Furthermore, the speed of the standard eigenaxis 
maneuver can be improved by 36% by delivering the optimal torque inputs directly to the 
reaction wheels themselves. These two software solutions reduce the reaction wheel 
advantage angle by 15% and 18%, respectively. The CMG performance gap is closed and 
surpassed in many cases. 

Due to their high torque capacity, CMGs dominate reaction wheels for small 
angle maneuvers. This is because a spacecraft using CMGs can accelerate to maximum 
angular velocity very quickly. Although CMGs have a lot more torque potential, they are 
exceedingly momentum-limited. A spacecraft using reaction wheels, on the other hand, is 
slower to accelerate, but has a much greater maximum velocity. Considering like-sized 
wheels, i.e., the same momentum capacity, this thesis shows that reaction wheels catch up 
and surpass CMGs due to their higher rate limits, provided the maneuver is large enough. 
The ramifications are particularly large for smallsats and cubesats as they require very 
small MEDs. For small satellites, this thesis has shown that reaction wheels have the 
ability to outperform like-sized CMGs at relatively small angle maneuvers. As wheel and 
spacecraft sizes decrease proportionally, reaction wheels catch up to CMGs even more 
quickly. Optimal control further improves this performance. 
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B. FUTURE WORK 

This thesis did not include an exhaustive approach at validating the models of the 
analytical effective eigenaxis of each of the optimal control models. Only 100 random 
data points were used for a single combination of spacecraft/wheel sizing. The results of 
this work should be validated across the entire space of reaction wheel and spacecraft 
sizes. Gimbal momentum/torque and gyroscopic coupling relationships should be 
accounted for into this “rule of thumb.” For instance, as body rate increases, the 
gyroscopic coupling increases. Most likely, the optimal maneuvers will use this as an 
advantage to widen the performance gap over standard maneuvers. 

The reaction wheel and CMG data should be analyzed further to account for 
variations in momentum size and torque output. To further understand the CMG 
momentum and torque relationship, research must be conducted to determine why a 
1:1 relationship between momentum and output torque is so commonly used. 

The different control approaches should also be explored from non-rest initial and 
final conditions. Given some calculable motion within a system, the optimal methods 
should be able to capitalize on this motion for further gains over the heritage maneuvers. 

Optimal control of the CMG should be explored. The optimal maneuver, by its 
very nature, should be able to avoid the singularities altogether. This would allow the 
unused capacity of the CMG momentum space to be accessed without complicated 
singularity-avoidance algorithms. One consideration during this expansion would be the 
confidence level in the spacecraft inertia and states. If knowledge of these uncertainties is 
poor, then the optimal control approach would still have a pseudo-inverse limited and 
singularity-ridden feedback system to deal with. Therefore, methods such as those 
described in [31] could be tried. 

Further analysis should be done on the momentum and torque envelope 
visualization method for optimal control maneuvers. Looking at these envelopes closely, 
it becomes obvious that the torque and momentum vectors tend toward the vertices and 
edges of the envelopes. If this space is studied in detail, it should be possible to predict 
the optimal slew times much more accurately. Just as the optimal eigenaxis-restricted 
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maneuver is a bang-eoast-bang optimal path, eertain eharaeteristies of optimal maneuvers 
eould be derived for maneuver size and eigenaxis orientation. This eould allow quiek on¬ 
board proeessing of semi-optimal maneuvers, without aetually eomputing the exact 
computer-intensive optimal path. A controller much like the eigenaxis quaternion 
feedback could potentially be developed to drive the system toward the appropriate 
vertices and edges of the momentum and torque envelopes. 

A comparison should be done with various reaction wheel and CMG 
arrangements. For example, if it is known that the pseudo-inverse limited time-optimal 
control logic is going to be used, then it may be possible to adjust the skew angle to 
maximize its envelope volume, thereby creating a larger momentum/torque radius (per 
Chapter V, Section G). Additionally, sizing wheels with different arrangements of wheel 
sizes could be considered. For example, a skew arrangement of five, six, or possibly ten 
smaller wheels may be more capable than four larger wheels, not to mention more 
redundancy. The shape would have more flat facets, but less loss. Like a cube and a 
regular octahedron of similar volumes, the octahedron allows a larger inscribed sphere. 
Mixed configurations of larger and smaller wheels could also be considered. Finally, 
hybrid reaction wheel and CMG configurations could be considered to maximize torque, 
yet maintain a high momentum envelope. 

Explore the impact of external disturbances on optimal maneuvers. Integrate 
gravity gradient, solar pressure, and aerodynamic torques into the optimal control 
problem. Compare with the eigenaxis maneuvers. Most likely, the optimal maneuvers 
will be able to exploit the external disturbances to widen the performance gap between 
the shaped eigenaxis trajectory and the time-optimal trajectories. 
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APPENDIX A. OPTIMAL CONTROL PROBLEM IN DESIGNER 
UNITS 
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APPENDIX B. MARKET RESEARCH 


A. REACTION WHEELS 


Manufacturer 

Model 

Mass 

(kg) 

1 ^ 

*^max 

(N-m-s) 

T 

^max 

(N-m) 

Reference 

Bradford Engineering 

W18 

4.95 

18 

0.248 

[ 32 ] 

Bradford Engineering 

WISE 

5.5 

22 

0.248 

[ 32 ] 

Bradford Engineering 

WISES 

5.85 

25 

0.248 

[ 32 ] 

Bradford Engineering 

W45 

6.45 

40 

0.248 

[ 32 ] 

Bradford Engineering 

W45HT 

6.5 

22 

0.403 

[ 32 ] 

Bradford Engineering 

W45E 

6.65 

55 

0.248 

[ 32 ] 

Bradford Engineering 

W45ES 

6.86 

70 

0.248 

[ 32 ] 

Dynacon MSCI 

Micro Wheel 4000 

5.2 

4 

0.05 

[ 33 ] 

Dynacon MSCI 

Micro Wheel 1000 

1.4 

1.1 

0.03 

[ 34 ] 

Dynacon MSCI 

Micro Wheel 200 

1 

0.18 

0.03 

[ 35 ] 

Honeywell 

HR0610 

5 

12 

0.055 

[ 36 ] 

Honeywell 

HR0610 

3.6 

4 

0.055 

[ 36 ] 

Honeywell 

HR12-12 

6 

12 

0.1 

[ 37 ] 

Honeywell 

HR 12-25 

7 

25 

0.15 

[ 37 ] 

Honeywell 

HR12-50 

9.5 

50 

0.2 

[ 37 ] 

Honeywell 

HR14-25 

7.5 

25 

0.1 

[ 37 ] 

Honeywell 

HR14-50 

8.5 

50 

0.15 

[ 37 ] 

Honeywell 

HR14-75 

10.6 

75 

0.2 

[ 37 ] 

Honeywell 

HR 16-50 

6 

50 

0.1 

[ 37 ] 

Honeywell 

HR16-75 

7 

75 

0.15 

[ 37 ] 

Honeywell 

HR 16-100 

9.5 

100 

0.2 

[ 37 ] 

E3 

RWA-15 

14 

20 

0.68 

[ 38 ] 

E3 

MWA-50 

10.5 

67.8 

0.07 

[ 39 ] 

Maryland Aerospace 

MAI-101 

0.64 

0.0011 

0.000635 

[ 40 ] 

Maryland Aerospace 

MAI-201 

0.73 

0.0108 

0.000625 

[ 41 ] 

Maryland Aerospace 

MAI-300 

0.317 

0.0076 

0.000625 

[ 42 ] 

NASAGSEC 

SMEX RW 

3.3 

4 

0.14 

[ 43 ] 


Table 7. Market Research Reaction Wheels 
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Manufacturer 

Model 

Mass 

(kg) 

h-max 

(N-m-s) 

T 

^max 

(N-m) 

Reference 

Rockwell Collins 

RSI 01-5/15 

0.6 

0.04 

0.005 

[ 44 ] 

Rockwell Collins 

RSI01-5/28i 

0.7 

0.12 

0.005 

[ 44 ] 

Rockwell Collins 

RSI 12-75/60 

4.85 

12 

0.075 

[ 44 ] 

Rockwell Collins 

RSI 4-75/60 

3.7 

4 

0.075 

[ 44 ] 

Rockwell Collins 

RSI 15-215/20 

7.7 

15 

0.215 

[ 44 ] 

Rockwell Collins 

RSI 45-75/60 

7.7 

45 

0.075 

[ 44 ] 

Rockwell Collins 

RSI 12-220/45 

6 

12 

0.22 

[ 44 ] 

Rockwell Collins 

RSI 18-220/45 

6.45 

18 

0.22 

[ 44 ] 

Rockwell Collins 

RSI 25-220/45 

7.15 

25 

0.22 

[ 44 ] 

Rockwell Collins 

RSI 50-220/45 

9.2 

50 

0.22 

[ 44 ] 

Rockwell Collins 

RSI 30-280/30 

9.2 

30 

0.28 

[ 44 ] 

Rockwell Collins 

RSI 68-170/60 

9.5 

68 

0.17 

[ 44 ] 

Rockwell Collins 

RSI 68-170/60 

9.5 

68 

0.17 

[ 44 ] 

Rockwell Collins 

RDR 23-0 

5.5 

23 

0.09 

[ 44 ] 

Rockwell Collins 

RDR 57-0 

7.6 

57 

0.09 

[ 44 ] 

Rockwell Collins 

RDR 68-3 

7.6 

68 

0.075 

[ 44 ] 

Rockwell Collins 

MWI 100-100/100 

16.5 

100 

0.1 

[ 44 ] 

Rockwell Collins 

MWI 30-400/37 

15.3 

30 

0.4 

[ 44 ] 

Sinclair 

Picosatellite RW 

0.09 

0.007 

0.001 

[ 45 ] 

Sinclair 

Picosatellite RW 

0.12 

0.01 

0.001 

[ 45 ] 

Sinclair 

Nanosatellite RW 

0.185 

0.03 

0.002 

[ 45 ] 

Sinclair 

Microsatellite RW 

0.225 

0.06 

0.005 

[ 45 ] 

Sinclair 

Microsatellite RW 

0.225 

0.06 

0.025 

[ 45 ] 

Sinclair 

Microsatellite RW 

0.97 

1 

0.1 

[ 45 ] 

SunSpace 

SunStar 

1.986 

0.65 

0.05 

[ 46 ] 

Surrey Space SSTL 

Micro Wheels lOSP-M 

0.96 

0.42 

0.011 

[ 47 ] 

Surrey Space SSTL 

Micro Wheels lOOSP-0 

2.6 

1.5 

0.11 

[ 47 ] 

Surrey Space SSTL 

Smallsat Wheels 200SP-M/O 

5.2 

12 

0.24 

[ 48 ] 

Surrey Space SSTL 

Smallsat Wheels 200SP-0 

5.2 

12 

0.24 

[ 48 ] 

Teldix 

20 cm 

3.4 

6.5 

0.2 

[ 5,p.398] 

Teldix 

26 cm 

6 

20 

0.2 

[ 5,p.398] 

Teldix 

35 cm 

8 

80 

0.2 

[ 5,p.398] 

Teldix 

50 cm 

12 

300 

0.3 

[ 5,p.398] 

Teldix 

60 cm 

37 

1000 

0.6 

[ 5,p.398] 


Table 8. Market Researeh Reaction Wheels (cont.) 
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B. 


CMGS 


Manufacturer 

Model 

Mass 

(kg) 

1 - 

*^max 

(N-m-s) 

T 

max 

(N-m) 

■ 

^max 

(rad/s) 

Reference 

Astrium 

CMG 4-6S 

13 

4 

6 

1.5 

[ 49 ] 

Astrium 

CMG 15-45S 

18.4 

15 

45 

3 

[ 50 ] 

Astrium 

CMG 10-30 

18.4 

10 

30 

3 

[ 51 ] 

Goodrich/Ithaco 

Violet CMG 

1.28 

0.31 

0.3100 

1 

[ 52 ] 

Honeywell 

M50 

33.1 

75 

75 

1 

[ 53 ] 

Honeywell 

M95 

38.6 

129 

129 

1 

[ 53 ] 

Honeywell 

M160 

44 

217 

217 

1 

[ 53 ] 

Honeywell 

M225 

54 

305 

305 

1 

[ 53 ] 

Honeywell 

M325 

61.2 

441 

441 

1 

[ 53 ] 

Honeywell 

M325D 

61.2 

441 

441 

1 

[ 53 ] 

Honeywell 

M715 

89.8 

969 

969 

1 

[ 53 ] 

Honeywell 

M600 

81.6 

813 

813 

1 

[ 53 ] 

Honeywell 

Ml 300 

125 

2000 

2000 

1 

[ 53 ] 

Honeywell 

M1400 

132 

2000 

2000 

1 

[ 53 ] 

Surrey Space SSTL 

BILSAT-1 CMG 

2.2 

0.28 

0.0559 

0.20 

[ 54 ] 

Swampsat 

Swampsat 

0.5 

8.00E-04 

8.00E-04 

1 

[ 55 ] 


Table 9. Market Researeh Control Moment Gyroscopes 
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